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ABSTRACT 


This thesis describes the development and demonstrates the use of 
SUBTRAN, a_ submarine transient detection computer simulation model. 
SUBTRAN provides, at the microcomputer level, a framework that can be used 
to investigate how transient detection and false transient detection opportunities 
affect the expected time to steady-state (continuous) detection. Monte Carlo 
methods are employed to simulate a submarine versus submarine passive 
acoustic detection search scenario. The scenario terminates when steady-state 
detection occurs. Detection is modeled using a signal excess threshold crossing 
model. Random fluctuations in the acoustic signal excess are modeled using a 
Lambda-Sigma Jump process. Both submarines are assumed to have a fixed 
speed and are constrained within a defined search area. The "target" submarine 
is assumed to be unable to counter-detect the "searching" submarine. Transient 
signals and false transient signals are determined by independent Poisson 
processes. Summary statistics for the times to detection are provided as output 


of the simulation. 
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I. INTRODUCTION 


Passive detection of a submarine can be divided into two mutually exclusive 
events: steady-state detection and transient detection. Steady-state detection 
involves the passive acoustic monitoring of a submarine’s continuous operating 
sounds. With today’s submarine quieting capabilities, steady-state detection 
typically requires short ranges. On the other hand, transient detection is the 
detection of a submarine’s intermittent sounds that are produced during 
infrequent operations. Due to the high signal levels generated by many of these 
infrequent operations, transient detection can occur at very long ranges. 

Of great concern to an evading submarine is the ability to freely move about 
without its actions being monitored by a searching submarine. This monitoring is 
generally based on an evading submarine’s continuous operating sounds. 

Due to the typical short ranges involved, the chance of a_ searching 
submarine being able to monitor an evading submarine’s continuously generated 
signals is quite small when a large search area is involved. However, an evading 
submarine’s transient signals offer long-range detection opportunities. The 
exploitation of these opportunities may improve the ability of the searching 
submarine to gain steady-state detection. 

A drawback to the use of long range detection opportunities occurs in a false 
transient environment. In any search scenario, false detections can lead the 
searcher away from the actual target position. 

The purpose of this thesis is to develop a means with which to address the 
issue of whether or not transient detection opportunities can be utilized by a 


searching submarine in a false transient environment to increase its chance of 


gaining steady-state detection on a target submarine. This will be accomplished 


by use of a Monte Carlo simulation, SUBTRAN. 
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Il. BACKGROUND 


A. GENERAL DESCRIPTION OF SIMULATION MODEL 

SUBTRAN is a discrete event simulation of two submarine platforms in a 
area search scenario. One platform is designated the "target and one the 
"searcher". During each discrete event time step, the time to the next discrete 
event is computed, steady-state detection is assessed, target and searcher 
positions are updated and the simulation time is advanced. A simulation run is 
completed when steady-state detection occurs or the user-defined maximum 
search time is reached. Propagation loss curves and figure of merits are used to 
derive the deterministic element of signal excess. A Lambda-Sigma Jump process 
is used to describe stochastic variations in signal excess, and a threshold 
crossing model is used to evaluate total signal excess for the occurrence of 
detection. Target characteristics incorporated in the simulation model include: 
target’s course, target’s speed, and target’s transient occurrence rate. Searcher 
characteristics incorporated in the simulation model include the searcher’s course 
and speed. Target depth, relative to the searcher, is assumed constant. False 
transients are modeled in the same fashion as target transients by using a 


transient occurrence rate. 


B. SEARCH SCENARIO 

The search scenario modeled in SUBTRAN is a parallel sweep search of a 
square area whose boundaries are known by both the target and searcher. At 
the beginning of each simulation replication, the target is positioned such that 


its X component of position and Y component of position are independent 


uniform random variables on the interval (0,L). L is the length dimension of the 
search area in nautical miles. The searcher’s position, at the beginning of each 
replication, is in the lower left-hand corner of the search area. In order to 
maximize steady-state search effort, the searcher is positioned a distance equal 
to one-half a track width from both the bottom and left-hand search area 
boundaries. The search area is oriented such that the top area boundary is to 
the north for course reference. A typical search area, annotated with initial 


searcher and target positions, is displayed in Figure 2.1. 


x 


TARGET’S POSITION 


SEARCHER’S POSITION 





Figure 2.1. Typical Search Area. 


C. DETECTION MODEL 

The basic detection model used in SUBTRAN converts user-defined 
propagation loss curves, figure-of-merits (FOM), and Lambda-Sigma Jump 
process parameters into an assessment of detection occurrence by evaluating 
total signal excess using a threshold crossing detection model. Two types of 
detection are modeled: steady-state detection and transient detection. Due to the 
current trend in submarine quieting, convergence zone detection of the target’s 
steady-state sound signature is assumed not to be possible. Convergence zone 
detection of the target’s transient noise signature is assumed to be possible but 
only to the third convergence zone. 

1. Propagation Loss 

The propagation loss experienced by sound as it travels through the 
water is determined by many factors. Frequency of the sound, receiver depth, 
source depth, ocean bottom type, and sound velocity profile all have an effect on 
the intensity of the sound received from an underwater source. In the area 
search scenario modeled in SUBTRAN, most of these factors remain relatively 
static with the exception of the frequency of the sound. Since steady-state and 
transient sounds are produced differently, they are likely to have different 
frequency components. As a result, two propagation loss curves are modeled in 
SUBTRAN, a steady-state propagation loss curve and a transient propagation 
loss curve. 

SUBTRAN allows the user to input up to 50 values of range and 
corresponding propagation loss for each of the two types of sound. From this 
input, the model constructs propagation loss curves using straight line 
approximation between successive points. One curve is constructed for steady- 


state signal loss and one is constructed for transient signal loss. Since 


convergence zone detection of the target’s steady-state sound signature is 
assumed not to be possible, only the direct path range on the steady-state 
propagation loss curve is used. 
2. Figure-of-Merit 
FOM, defined for passive sonar as the maximum propagation loss which 
still results in a 50 percent probability of signal detection [Ref. l:p. 49], is 


modeled in SUBTRAN under the following assumptions and conditions: 


- The target cannot counter-detect the searcher. 
- The speeds of the target and searcher are held constant. 


- There is only one type of steady-state sound source emitted by the target 
and it is omni-directional. 


- There is only one type of transient sound source emitted by the target and 
it is omni-directional. 


- The ambient noise level is constant throughout the simulation. 


Under these assumptions and conditions, target FOM can _ be 
summarized by two constant values: steady-state FOM and transient FOM. 
SUBTRAN requires the user to determine the values of the steady-state and 
transient FOM for the target based on the search scenario being simulated. 

3. Signal Excess 

Signal excess (SE) is modeled in SUBTRAN as consisting of a 
deterministic part and a stochastic part. The deterministic part of SE is defined 
as the algebraic difference between FOM and propagation loss. The stochastic 
part of signal excess is used to describe the random fluctuations that occur in 
signal excess with respect to time. A Lambda-Sigma Jump process is used to 


model these random fluctuations. 


4. Lambda-Sigma Jump Process 

A Lambda-Sigma Jump process [Ref. 2:p. 8] is used to describe the 
random fluctuations in signal excess. This process models the time between 
fluctuations as well as the magnitude of the fluctuations. For the Lambda-Signal 
Jump process, the time between fluctuations is exponentially distributed with 
rate parameter lambda and the magnitude of the fluctuations is normally 
distributed with mean zero and standard deviation sigma. 

The user inputs values for the parameters lambda and_ sigma. 
Appropriate values for these parameters are subjective and must be based on the 
experience of the user. A guide in the estimation of lambda and sigma, based on 


various sources, 1s [Ref. 3 and Ref. 4]: 


- Lambda varies between 2 and 60 hours". 


- Sigma varies between 3 and 9 GB. 


Other sources containing information helpful in determining appropriate values 
for lambda and sigma are listed by Tehan [Ref. 4]. 
5. Threshold Crossing Model 
A threshold crossing model [Ref. 2:p. 7] is the final step in the overall 
effort to model detection. This model predicts that detection occurs when signal 
excess is equal to or exceeds some given threshold level. The most common 
threshold level, and the one employed in SUBTRAN, is zero. In other words, 
detection occurs whenever signal excess is greater than or equal to zero. 
6. Detection Model Employment. 
The detection model outlined in the previous sections is implemented in 
SUBTRAN by extracting the current detection range from the propagation loss 


curve and comparing it to the range between the target and the searcher. If the 


detection range is greater than or equal to the range between the target and the 
searcher, detection occurs. 

The current detection range is determined by considering an adjusted 
value of the FOM. This adjusted value is obtained by adding the user-defined 
value of FOM and the current signal excess fluctuation value. The range at 
which this adjusted FOM value and the propagation loss curve are equal is the 
current detection range. 

7. Steady-State Detection Assessment 

SUBTRAN is a discrete event simulation. A discrete event is defined as 
any event whose occurrence could change the ability of the searcher to gain 
steady-state detection on the target. The discrete events modeled in SUBTRAN 
are: target course change, target boundary reflection, searcher course change, 
target transient signal occurrence, false transient signal occurrence, and signal 
excess fluctuation occurrence. 

Since SUBTRAN is a discrete event simulation, steady-state detection 
must be assessed for the time period until the next discrete event. To do this, 
the current values of searcher and target course and position are used together 
with the user-defined speed values to define the positions of the searcher and 
target at any future time. These positions are then compared to find the time it 
would take the range between the searcher and target to decrease to a value 
less than or equal to the current steady-state detection range. If the time 
calculated is less than or equal to the time to the next discrete event, steady- 
state detection occurs. Appendix E contains the steady-state detection assessment 


calculations. 


8. Transient Detection Assessment 

Detection of target transients by the searcher is assessed in a different 
manner than steady-state detection. The target’s transient sound signature is 
intermittent and when produced, is assumed to occur instantaneously. As a 
result of these characteristics, transient detection is assessed in SUBTRAN by 
comparing searcher and target positions at the time of a transient occurrence. If 
the range is within the transient detection range, transient detection occurs. 

Transient detection range is extracted from the transient propagation 
loss curve using the transient FOM. Since transient detection can occur out to 
the third convergence zone, the transient detection range consists of four range 
brackets: a direct path range bracket and three convergence zone annulus range 


brackets. 


D. TARGET CHARACTERISTICS 
The target characteristics modeled in SUBTRAN include: speed, course, and 
transient emission rate. Target depth relative to the searcher is assumed to be 
constant. In order to model its worst case vulnerability to detection, the target 
does not have the ability to detect and thus avoid the searcher. 
1. Target Speed 
Target speed is held constant throughout the simulation at a value 
specified by the user. 
2. Target Course 
The search scenario modeled in SUBTRAN involves a randomly moving 
target which is constrained to a defined search area by boundary reflection. To 


incorporate differing target course and time to course change distributions, 


SUBTRAN uses the inverse probability integral transformation [Ref. 5:pp. 20- 
21a] to model target movement. 

The inverse probability integral transformation method involves 
selecting values of a desired random variable whose cumulative distribution 
function (CDF) is known by using random variables drawn from a uniform (0,1) 
distribution. This method is illustrated in Figure 2.2. A random variable drawn 
from a uniform (0,1) distribution is entered on the Y-axis of the CDF 
representing the random variable desired. The corresponding X-axis value is the 
associated random variable of the given distribution. 

SUBTRAN allows the user to input up to 50 values for both the target’s 
course CDF and the time to course change CDF. The corresponding CDF curves 
are constructed using straight-line approximation between successive points. 

Due to the randomness of the target's movement, a procedure exists to 
constrain the target to the search area. Using its current course and speed, 
the position of the target is evaluated for the period until its next course change. 
If during this time the target encounters the search area boundary, its course is 
reflected back into the search area. The angle of reflection is equal to the angle 
of incidence between the search area boundary and the target’s incoming course. 
An example of target boundary reflection is shown in Figure 2.3. 

3. Target Transients 

Target transient occurrence is modeled in SUBTRAN as a Poisson 
process with a user-defined rate parameter. Only one type of transient sound is 
assumed to be produced by the target. In addition, when the transient does 


occur, it is assumed to be instantaneous. 
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Figure 2.2 Inverse Probability Integral Transformation Method. 
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E. SEARCHER CHARACTERISTICS 
Searcher characteristics modeled in SUBTRAN include: speed and course. 
The depth of the searcher relative to the target is assumed constant. 
1. Searcher Speed 
Searcher speed is held constant throughout the simulation at a value 
specified by the user. 
2. Searcher Course 
Searcher movement, unlike target movement is deterministic. The 
movement modeled in SUBTRAN consists of a general search movement and a 
transient prosecution movement. 
The general search movement is an area search using parallel sweeps. 
The spacing between sweeps is defined as the track width, TW. TW can be user 
specified or defaulted. When defaulted, TW is set equal to twice the value of the 
mean steady-state detection range. To maximize the amount of search effort 
applied to the search area, the searcher is programmed to stay a distance of 
(0.5) x (TW) from the area boundary. This distance constraint forces an 
adjustment to the searcher’s movement as it approaches the end of the area 
search. Specifically, after each vertical leg of the search, the position of the 
search is checked. If the horizontal distance between the searcher and the search 
area boundary that lies ahead of the searcher is less than (1.5) x (TW), the next 
horizontal leg of the search is adjusted. The adjustment that occurs is a 
shortening of the distance the searcher normally travels. The normal distance of 
the horizontal leg of the search is TW. The adjusted distance, TW,,, is given by 


the following expression: 


TW,,, = TW x INT(L/TW) 
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where, 


L length dimension of the search area 


INT(L/TW) 


integer part of L/TW 


Figure 2.4 illustrates the track width adjustment and displays the 
resulting return area search that occurs. Maximum search area coverage 18 
achieved when L/TW is an integer value. 

The prosecution phase of searcher movement involves suspending or 
terminating the current search, investigating a transient detection, and returning 
to the general area search if an additional transient detection does not occur. 

Upon detection of a transient, the searcher either suspends the general 
area search or terminates a prosecution movement, depending on what search 
phase is being conducted, and commences a new prosecution movement. A 
prosecution movement involves steering the course on which the transient was 
detected until the search area boundary is approached. When the distance 
between the searcher and the search area boundary along the searcher’s track 
reduces to (0.5) x (TW), the searcher reverses direction and continues down the 
reciprocal bearing until he returns to the point from which the prosecution 
movement was commenced. From this point, the searcher returns to the position 
where the general area search was suspended and resumes the general area 
search. If at any time during the prosecution movement phase another transient 
is detected, the searcher terminates the current movement phase and commences 


a new one based on the most recent transient detection. 
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F. FALSE TRANSIENT MODEL 

False transients are included in SUBTRAN to represent transient sound 
sources that the searcher may mistakenly identify as target transients. The 
characteristics of false transients that are modeled in SUBTRAN include: the 
rate of occurrence, the position at time of occurrence, and the ability of the 
searcher to detect the occurrence. 

False transient occurrence is modeled as a Poisson process. As a result, the 
rate of false transient occurrence is exponentially distributed with a user-defined 
rate parameter. 

In an effort to best represent natural conditions, the position from which a 
false transient is emitted is assumed to be uniformly distributed over the search 
area. 

The ability of the searcher to detect the occurrence of a false transient is 
assessed in the same manner as the detection of real transients. Specifically, the 
range between the searcher and the position that the false transient was emitted 
from is compared to the detection ranges extracted from the transient 
propagation loss curve. The FOM for false transients is assumed to be the same 


as the FOM for real transients. 


G. SUMMARY OF IMPORTANT MODEL ASSUMPTIONS 
The major assumptions made in the development of SUBTRAN include: 
- The search area is square and its boundaries are known to both the target 
and searcher. 
- No convergence zone steady-state detection. 
- Third convergence zone transient detection. 
- Target cannot counter-detect the searcher. 
- Only one type of target transient sound. 
15 


Target transient sounds are instantaneous. 
FOM is constant. 


Random fluctuations in signal excess are described by a Lambda-Sigma 
Jump process having the values of lambda and sigma specified by the user. 


Target and false transients are independent Poisson processes. 
False transients sounds have the same FOM as target transient sounds. 


Target position is constrained to the search area by reflection. 
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HI. DESCRIPTION OF SUBTRAN 


A. PROGRAM CODE AND STRUCTURE 

SUBTRAN is coded in Fortran 77 [Ref. 6] and compiled using IBM 
Professional Fortran. It is designed to be implemented on a personal computer. 
The program is structured in modular form with major tasks such as boundary 
reflection, assessment of steady-state detection, and searcher movement being 
accomplished by separate subroutines. This type of structuring facilitates easy 
alteration of major program functions by a knowledgeable programmer. However, 
a compiled version of the program is all that is needed by most users. The 
compiled version was designed to allow a user, who accepts the assumptions 
made in this thesis, to change a select set of input parameters in an input data 


file and not have to perform any programming. 


B. MAIN PROGRAM AND ASSOCIATED SUBROUTINES 

The programming code for SUBTRAN is broken up into 19 parts: a main 
program and 18 supporting subroutines. The interrelationship of the various 
parts is shown in Figure 3.1. 

1. Main Program 

The main program, SUBTRAN.FOR, controls the call of the various 

subroutines. Specific functions performed in SUBTRAN.FOR include: 
initialization of counters, initial positioning of the searcher and _ target, 
calculation of the time to the next discrete event, control of the simulation time 
clock, storage of the times to steady-state detection, and control of the number of 


simulation replications. 
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Figure 3.1. Main Program and Subroutine Interrelation. 
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2. Input Subroutine 
The input subroutine, IN.FOR, is responsible for reading the input data 
file and passing the input parameters to the main program. Other functions 
performed in IN.FOR include: printing the input parameters to the user's 
terminal for verification and printing the input parameters to the output data 
file SUBTRAN.OUT. 
3. Propagation Loss Subroutine 
The propagation loss subroutine, PLOSS.FOR, is_ responsible for 
calculating the steady-state and transient detection ranges. Major parameters 
passed to PLOSS.FOR by the main program include: the propagation loss curve 
values and the adjusted FOM. 
4. Target Course Subroutine 
The target course subroutine, TGCR.FOR, is responsible for calculating 
the target’s next course and time to course change. Major parameters passed to 
TGCR.FOR by the main program include: the target’s course and time to course 
change CDF curve values and a random number seed. 
5. Target Transient Subroutine 
The target transient subroutine, TGTR.FOR, is_ responsible for 
calculating the time to the target’s next transient. Major parameters passed to 
TGTR.FOR by the main program include: the target’s transient occurrence rate 
parameter and a random number seed. 
6. Searcher Course Subroutine 
The searcher course subroutine, SRCR.FOR, is responsible for 
calculating the searcher’s next course and time to course change during general 
area search movement. Major parameters passed to SRCR.FOR by the main 


program include: the searcher’s position and speed, track width, length 


fg 


dimension of the search area, and searcher’s course flag. The searcher’s course 
flag indicates which part of the general area search is to be performed next. 
7. Transient Course Subroutine 
The transient course subroutine, TRCR.FOR, is responsible for 
calculating the searcher’s next course and time to course change during transient 
prosecution movement. Major parameters passed to TRCR.FOR by the main 
program include: the searcher’s position and speed, track width, length 
dimension of the search area, and searcher’s course flag and transient course 
flag. The searcher’s transient course flag indicates two things: whether the 
searcher is performing a general area search movement or transient prosecution 
movement, and if a transient prosecution movement is being performed, which 
leg of the movement is to be performed next. 
8. Steady-State Detection Subroutine 
The steady-state detection subroutine, SSDET.FOR, is responsible for 
calculating the time to Syed ane detection. Major ae passed to 
SSDET.FOR by the main program include: the position of the searcher and 
target, the course and speed of the searcher and target, and the steady-state 
detection range. 
9. Transient Detection Subroutine 
The transient detection subroutine, TRDET.FOR, is responsible for 
assessing the occurrence of target transient and false transient detection. Major 
parameters passed to TRDET.FOR by the main program include: the position of 
the searcher and the transient emission, and the transient detection ranges. 
10. Lambda-Sigma Jump Subroutine 
The Lambda-Sigma Jump subroutine, LSJ.FOR, is responsible for 


calculating the magnitude of the next signal excess fluctuation as well as the 
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time to occurrence. Major parameters passed to LSJ.FOR by the main program 
include: the Lambda-Sigma Jump error process rate parameter and standard 
deviation, and a random number seed. 
11. False Transient Subroutine 
The false transient subroutine, FALTR.FOR, is_ responsible for 
calculating the position of and time to the next false transient. Major parameters 
passed to FALTR.FOR by the main program include: the false transient rate 
parameter, the length dimension of the search area, and a random number seed. 
12. Position Subroutine 
The position subroutine, POSIT.FOR, is responsible for updating the X 
and Y positions of the searcher and target. Major parameters passed to 
POSIT.FOR by the main program include: the previous X and Y positions of the 
searcher and target, the course and speed of the searcher and target, and the 
time passed since the last update. 
13. Boundary Reflection Subroutine 
The boundary reflection subroutine, REFLCT.FOR, is responsible for 
calculating the targets next boundary reflection course and the time to reflection. 
Major parameters passed to REFLCT.FOR by the main program include: the 
target’s position, the target’s course and speed, and the length dimension of the 
search area. REFLECT.FOR is used by the TRACK.FOR subroutine in addition 
to the main program. 
14. Track Subroutine 
The track subroutine, TRACK.FOR, is responsible for calculating the 
course and time necessary for the searcher to reach a given point in the search 
area. Major parameters passed to TRACK.FOR by the main program include: the 


positions of the searcher and the point to be reached, the speed of the searcher, 
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and the length dimension of the search area. TRACK.FOR is used by the 
TRCR.FOR subroutine in addition to the main program. 
15. Output Subroutine 

The output subroutine, OUT.FOR, is responsible for writing the output 

data to the user’s terminal and to the two output data files. 
16. Miscellaneous Subroutines 

Four miscellaneous subroutines are called by the main program and 
other subroutines: HSCRN.FOR, LRNDPC.FOR, LNORPC.FOR, and 
LGAMPC.FOR. 

The subroutine HSCRN.FOR is used by the main program, the input 
subroutine, and the output subroutine to stop the output display on the user’s 
terminal. 

The subroutines LRNDPC.FOR, LNORPC.FOR, and LGAMPC.FOR are 
used by the main program and various subroutines to generate random numbers. 
All three are an Me version of the subroutines, with the same names, 
developed by Lewis, Orav, and Uribe [Ref. 7]. 

LRNDPC.FOR is used to generate uniform (0,1) random variables, 
LNORPC.FOR is used to generate normal (0,1) random variables, and 


LGAMPC.FOR is used to generate exponential (1) random variables. 


C. PROGRAM OUTPUT 

The output provided by SUBTRAN consists of: a_ verification of input 
parameters, various output statistics, and the specific times to detection. The 
output is written to two separate output files named SUBTRAN.OUT and 


TIME.OUDT. 
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1. Verification of Input Parameters 

Ail parameters entered via the input data file are printed in the output 
file SUBTRAN.OUT. This allows user verification of input as well as a hardcopy 
display. 

2. Various Output Statistics 

SUBTRAN performs elementary statistical analysis on the times to 
detection and provides the results in the output data file SUBTRAN.OUT. The 
statistics include: mean and standard deviation of the times to detection, 95% 
confidence bound for the mean, estimate of the mean time to detection assuming 
the times to detection are exponentially distributed, R’ confidence factur for the 
estimated mean time to detection, and the cumulative distribution function of 
the times to detection in tabular form. Calculations are performed on the times 
to detection conditioned on no detection at time zero. In other words, runs with 
detection time zero were disregarded. The number of simulation runs with a zero 
detection time is included in the output summary. 

The mean and standard deviation are computed using normal methods. 
For the replications in which steady-state detection does not occur, T,,,, 18 used 
as the detection time for the calculation of the mean and standard deviation. 
This results in a lower bound on the actual mean time to detection and standard 
deviation. The number of times to detection that are set to T,,,, 18 provided in 
the output summary. 

The 95% confidence bound for the mean is calculated using: 


95% bound = o/N” 


where, 
o = standard deviation of the times to detection. 


N = number of simulation replications. 
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The estimated mean time to detection, under the assumption that the 
times to detection are exponentially distributed, is calculated using straight line 
regression of the natural log of the empirical survival function on the ordered 
times to detection [Ref. 8:pp. 604-607]. The survival function of an exponential 
distribution when plotted on a natural log scale is a straight line. The 
magnitude of the slope of the straight line is equal to the rate parameter of the 
exponential distribution. The reciprocal of the rate parameter is the estimated 
mean time to detection. Replications in which steady-state detection did not 
occur are disregarded in the regression under the assumption that if the 
distribution is truly exponential, their inclusion would not affect the estimate of 
the slope. 

The R’ confidence factor for the estimated mean time to detection 
describes the proportion of the variation in the sample values of the natural log 
of the survival function that the regression estimate accounts for. R’ was 
calculated using normal methods [Ref. 8:p. 644]. ; 

The cumulative distribution function of the times to detection is 
empirically derived from the ordered times to detection [Ref. 9:pp. 11-16] and is 
provided in tabular form. An example of the output statistics provided by 
SUBTRAN in contained in Appendix D. 

3. Specific Times to Detection 

The specific times to detection are provided in the output data file 

named TIME.OUT. The times are ordered from shortest to longest with zero 


times removed and T,,,, value included as necessary. 
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IV. VALIDATION OF SUBTRAN 


Validation of a simulation model can be broken up into two major parts: 
validation that the computer code is functioning as expected and validation that 
the simulation model is an accurate assessment of the real world situation being 
modeled. 

With regard to SUBTRAN, validation that the computer code is functioning 
as expected was conducted using two separate techniques. First, the simulation 
was conducted with a fixed target and the time to detection was compared to a 
uniform distribution. Second, the simulation was conducted with a fixed searcher 
and the time to detection was compared to model results for a random tour 
target developed by Abd El-Fadeel [Ref. 10]. The results of these two tests were 
used to validate the proper simulation coding of searcher movement, target 
positioning, target movement, target boundary reflection, and steady-state 
detection assessment. 

Areas of program code not rigorously validated include searcher transient 
prosecution, transient occurrence, and false transient occurrence. In the case of 
transient occurrence and false transient occurrence, the program code is 
relatively straightforward and rigorous validation is not necessary. The program 
code for searcher transient prosecution is more complex and validation was 
conducted by reviewing searcher position data after simulation runs in which 
transients were allowed. The data reviewed indicated the searcher was 
performing transient prosecutions as modeled. 

Validation, in a strict sense, that the simulation model is an accurate 


assessment of the real world situation being modeled is not possible. There 
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exists little to no data that could be used in a rigorous validation of SUBTRAN. 
However, the results produced by the example runs of SUBTRAN lend credibility 


to the modeling technique used. The results clearly indicated intuitive results. 


A. UNIFORM DISTRIBUTION APPROXIMATION 

By altering the search scenario modeled in SUBTRAN, a_ uniform 
distribution of the times to detection should be achieved. This contention is 
obvious from the fact that at the start of each simulation replication the target 
is uniformly positioned in the search area. With the target’s speed set at zero, 
the time it takes the searcher to detect the target on one complete search of the 
area should be uniformly distributed on the open interval (0,T). Where T is the 
time it takes the searcher to make one complete search of the area. 

This uniform distribution approximation was tested on SUBTRAN using four 
different sets of input. The input values are summarized in Table IV.1. Figures 
4.1 through 4.4 graphically compare the estimated probability of non-detection as 
a function of search time with the theoretical probability of non-detection for the 
four cases. The theoretical probability of non-detection is indicated as a dotted 
line. The estimated mean time to detection (MTTD), its 95% confidence interval, 


and the theoretical MTTD are also included. 


TABLE IV.1. UNIFORM APPROXIMATION INPUT PARAMETERS. 


* 1000 replications performed in all cases. 





26 


The results of the uniform approximation test strongly indicate proper 
functioning of the simulation with regard to searcher movement and target 
positioning. One note worth mentioning deals with the slight graphical deviation 
in theoretical and estimated probability of non-detection shown in Figure 4.2, 
case 4. The graphical results indicate that detection in the simulated case occurs 
sooner than detection in the theoretical case at large values of time. The reason 
for this is explained by the geometry of the particular search. In the first three 
test cases the length dimension of the search area was evenly divisible by the 
searcher’s track width. As a result, a complete search of the area required no 
overlap of search effort. However, in the fourth case tested, the division of the 
length dimension by the track width did not result in a integer and overlap 
occurred on the last leg. This overlap caused the deviation in the CDF. Target 
detections that should have theoretically taken place on the last search leg 


occurred one leg earlier. 


B. RANDOM TOUR MODEL APPROXIMATION 

Shifting to a moving target and a stationary searcher, the results of 
SUBTRAN can be compared directly with the random tour model [Ref. 10:pp. 59- 
62]. The random tour model was developed using Monte Carlo simulation 
techniques and applies to the detection of a randomly moving target constrained 
to a square search area by boundary reflection. The model also assumes that the 
time between target course changes is exponentially distributed with rate 
parameter t and that the searcher is fixed in the center of the search area with 
a cookie-cutter detection radius. The output of the model estimates the 
probability of non-detection with respect to time for given values of searcher 
detection range (R), target speed (V), search area size (A), and the target course 
change rate parameter (t). The equation for the probability of non-detection 
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Figure 4.1. Uniform Approximation Case 1 - Probability of 
Non-Detection. 
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Figure 4.2. Uniform Approximation Case 2 - Probability of 
Non-Detection. 
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Figure 4.3. Uniform Approximation Case 3 - Probability of 
Non-Detection. 
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Figure 4.4. Uniform Approximation Case 4 - Probability of 
Non-Detection. 
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(PND), given that detection did not occur at time zero, proposed by the random 


tour model is: 
PND(t) = exp™ 
where, 
a = [(24.7RV’) / WA - 2R’)"*] x [1 - exp(-0.0841A"’/V)] 


Four SUBTRAN simulation runs were compared to the random tour model. 
The input parameters are summarized in Table IV.2. Figures 4.5 through 4.8 
graphically compare the estimated probability of non-detection as a function of 
search time with the random tour model estimate for probability of non-detection 
for the four cases. The random tour model estimate is indicated as a dotted line. 
The estimated MTTD, its 95% confidence interval, and the random tour model 
(RTM) estimate of MTTD are also included. 


TABLE IV.2. RANDOM TOUR MODEL APPROXIMATION INPUT 
PARAMETERS. 


* 1000 replications were performed in all cases. 





This test shows the close comparison between the results produced by 
SUBTRAN and those predicted by the random tour model. In all cases tested, 


SUBTRAN produced a slightly shorter time to detection than the random tour 
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Figure 4.5. Random Tour Model Approximation Case 1 - Probability of 
Non-Detection. 
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Figure 4.6. Random Tour Model Approximation Case 2 - Probability of 
Non-Detection. 
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Figure 4.7. Random Tour Model Approximation Case 3 - Probability of 
Non-Detection. 


10 nm 
100 nm 
10 nm/hr 


Estimated MTTD = 54+3 hrs 
RTM MTTD = 59 hrs 


= 
© 
= 
O 
Lu) 
to 
T 
= 
© 
= 
L. 
.@) 
a 
= 
- 
O 
aa 
os 


&0 
SEARCH TIME (HRS) 





Figure 4.8. Random Tour Model Approximation Case 4 - Probability of 
Non-Detection. 
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model. These results strongly validate the modeling of target movement, target 


boundary reflection, and computation of the occurrence of steady-state detection. 
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V. EMPLOYMENT OF SUBTRAN 


A. INTRODUCTION 

This chapter is designed to demonstrate the employment of SUBTRAN. Five 
separate search scenarios involving plausible input parameters were simulated 
using various combinations of target and false transient occurrence rates. The 
output statistics of these five test cases are presented in graphical form to 
emphasize the relationship between mean time to detection (MTTD), mean time 
between target transients (MTBTT), and mean time between false transients 
(MTBFT). In addition, the distributions of the times to detection are analyzed to 
assess the assumption that they belong to the exponential class of distributions. 

Two additional search scenarios, Excursion 1 and Excursion 2, were 


simulated to provide a further understanding of how false transient occurrences 


7 


affect MTTD. 
The numeric output statistics generated by SUBTRAN were used to develop 


the graphs displayed in this chapter. 


B. INPUT PARAMETERS 
The input parameters that remained common among all five test case and 


the two excursions are: 


- Steady-state sound propagation loss curve displayed in Figure 5.1. 


- Transient sound propagation loss curve equal to the steady-state sound 
propagation loss curve displayed in Figure 5.1. 


- Transient sound figure of merit set at 90 dB. 


- Target’s course cumulative distribution curve displayed in Figure 5.2. 
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- Target's time to course change cumulative distribution curve displayed in 
Figure 5.3. 


- Lambda-Sigma Jump model parameters. Lambda set at 3 hour’; Sigma at 6 
dB. 


- Track width flag set at 0. 


- Number of simulation replications set at 1000. 


The input parameters that varied over the five test cases and two 
excursions were: the mean steady-state detection range (R), the length dimension 
of the search area (L), the target’s speed (S,), the searcher’s speed (S,), and the 
maximum search time (T,,,,). The values of these parameters are summarized in 
Table V.1. 

The target and false transient occurrence rates, converted to MTBTT and 


MTBFT respectively, used in each test case and excursion are: 


- MTBTT’s of 1000, 200, 100, 75, 50, 10, 5 (in hours). 


- MTBFT’s of 400, 75, 10, 5, 1 (in hours). 


In addition, simulation runs were made in which no target and/or false 


transient sounds were allowed. 


C. OUTPUT STATISTICS AND ANALYSIS 

The MTTDs and associated 95% confidence intervals for the five test cases 
and two excursions are tabulated in Tables V.2 through V.8. 

Three dimensional graphical displays of the relationship between MTTD, 
MTBTT, and MTBFT for the five test cases and two excursions are contained in 
Figures 5.4 through 5.8, 5.14, 5.16. 

Graphical analyses of the probability of non-detection (PND) for the five test 


cases and two excursions are displayed in Figures 5.9 through 5.13, 5.15, 5.17. 
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Figure 5.1. Steady-State and Transient Sound Propagation Loss Curve. 
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Figure 6.2. Target’s Course Cumulative Distribution Curve. 
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Figure 5.3. Target’s Time to Course Change Cumulative Distribution 
Curve. 
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TABLE V.1. EMPLOYMENT INPUT PARAMETERS. 


Test Case 


* EX 2 involves deceptive false transient sounds. 





1. Test Cases 1 Through 5 
Test cases 1 through 5 involve variations in R, L, 8,, S,,and T,,,,. The 
variations of R, L, S,, and S, were chosen to encompass plausible search 
scenarios. T,,, was chosen, based on the search scenario involved, to minimize 
the number of simulation runs in which steady-state detection did not occur. 
a. Relationship Between MTTD, MTBTT, and MTBFT 
Figures 5.4 through 5.8 indicate the following three points: 
- eat decreases as MTBTT decreases throughout the range of MTBTT’s 
ested. 


- The "decrease’ trend is present in MTTD, as a function of MTBTT, for all 
values of MTBFT tested. 


- With respect to MTBFT, MTTD is relatively constant for a given MTBTT 
unless a high false transient occurrence rate is present (i. e. MTBFT < 5 
hours). 

These three points highlight the general results indicated by Test 


Cases 1 through 5. Specifically, target transient sounds occurring at relatively 


moderate rates aid the searching submarine in the steady-state search effort, 
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while false transient sounds only hamper the search effort when their occurrence 
rate is high. 

Figures 5.9 through 5.13 further emphasize the effects of target and 
false transient sounds on the searcher’s steady-state search effort. Three lines 
are plotted on each graph representing the PND for the three scenarios: no 
target or false transient sounds present, only target transient sounds present 
with a MTBTT of 5 hours, and both target and false transient sounds present 
with a MTBTT of 5 hours and a MTBFT of 1 hour. The improvement in the 
time to steady-state detection is indicated by comparing the line associated with 
no transient sounds and the line associated with only target transient sounds. 
The degradation in the time to steady-state detection caused by the introduction 
of false transients is indicated by comparing the line associated with only target 
transient sounds and the line with both target and false transient sounds. In all 
of the test cases the magnitude of the decrease in the time to steady-state 
detection realized by the introduction of target transient sounds with a MTBTT 
of 5 hours is greater than the magnitude of the increase in the time to steady- 
state detection realized by adding false transient sounds with a MTBFT of 1 
hour. 

b. Distribution of the Times to Detection 

Figures 5.9 through 5.13 can also be used to graphically assess the 
exponential ‘fit’ of the estimated PNDs. Theoretical PND curves, assuming the 
times to detection are exponentially distributed with a mean equal to the 
estimated mean time to detection, are plotted as a dashed lines for comparison. 
It is apparent that for the five test cases, PND can be closely approximated by 


an exponential distribution with a mean equal to the mean time to detection. 
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Additional support, in favor of the exponential approximation of the 
PND, was provided by examining the output statistics generated by SUBTRAN. 
The estimated mean of the times to detection and its associated 95% confidence 
interval, the estimated standard deviation of the times to detection, and the 
mean time to detection calculated using linear regression of the times to 
detection (see Chapter III.C.2 for explanation) were numerically compared for 
agreement. If the underlying distribution is truly exponential, all three values 
should be approximately equal. Of the 215 simulations involved in the five test 
cases, 205 simulations resulted in the difference between the estimated mean of 
the times to detection and the estimated standard deviation to be within the 
95% confidence interval associated with the mean, while 206 simulations 
indicated similar agreement between the estimated mean of the times to 
detection and the mean time to detection calculated using linear regression. The 
R’ confidence factor in all simulations was greater than 0.98 with the majority of 
the values greater than or equal to 0.99. 
2. Excursions 1 and 2 
In an effort to better understand when false transients play a major 
role in the steady-state search effort, two excursion scenarios were simulated. 
The first excursion involved a slow, almost stationary (speed of 0.25 kts/hr), 
target and a searcher with a large mean steady-state detection range. The 
second excursion involved ‘deceptive’ false transient sounds. 
In the first excursion, an attempt was made to test a search scenario in 
which false transient sounds should noticeably degrade the search effort. As 


presented in Chapter IV, the times to detection are uniformly distributed in a 
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TABLE V.2. TEST CASE 1 
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Figure 5.5. Test Case 2- Estimated MTTD, MYIBTT, and MTBFT 
Relationship. 


TABLE V.3. TEST CASE 2 - ESTIMATED MTTD VALUES WITH 95% CI. 
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Figure 5.6. Test Case 3 - Estimated MTTD, MTBTT, and MTBFT 
Relationship. 


TABLE V.4. TEST CASE 3 - ESTIMATED MTTD VALUES WITH 95% CI. 
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Figure 5.7. Test Case 4 - 


TABLE V.5. TEST CASE 4 - ESTIMATED MTTD VALUES WITH 95% CI. 
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Estimated MTTD, MTBTT, and MTBFT 
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Figure 5.8. Test Case 5 - 
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TABLE V.6. TEST CASE 5 - ESTIMATED MTTD VALUES WITH 95% CI. 
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Figure 5.9. Test Case 1 - Probability of Non-Detection. 
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Figure 5.10. Test Case 2 - Probability of Non-Detection. 
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Figure 5.11. Test Case 3 - Probability of Non-Detection. 
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Figure 5.12. Test Case 4 - Probability of Non-Detection. 
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Figure 5.13. Test Case 5 - Probability of Non-Detection. 


scenario involving a stationary target, a systematically moving searcher, and no 
transient sounds. However, by introducing target and false transient sounds to 
the scenario, the searcher is forced away from his systematic search by transient 
occurrence investigations. In essence, the search effort is degraded to a random 
search of the area. 

Comparing the mean time to detection in the uniform distribution case 


to the mean time to detection in the random search case we have: 


MTTDyarormm = A/ 4RV- MTTDeisson = A/ 2B 


where, 


area of search 


= radius of detection 


| 


effective search speed 
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Since the mean time to detection in the uniform distribution case is one 
half that in the random search case, the degradation in MTTD experienced with 
the introduction of false transient sounds should be greater than the degradation 
that occurred in the five test cases. Figure 5.14 indicates this contention. In 
addition, Figure 5.15, indicates the transition from a uniform search to a random 
search. 

In Figure 5.15, the curve associated with no false transient sounds has 
deviated from the exponential fit as is indicated by the bowing of the curve on a 
natural log scale. Finally, Figure 5.15 also refutes the trend, established in the 
five test cases, that the gains achieved by target transient sounds were greater 
than the losses experienced from the introduction of false transient sounds at 
nearly equal rates of occurrence. The curve associated with both target and false 
transient sounds shifts back to the position of the curve associated with no 
transients and at higher time values surpasses it. 

Excursion 2 deals with "deceptive' false transients. The main program of 
SUBTRAN was altered slightly to achieve a method for drawing the searcher 
away from the target in the event of a false transient sound emission. In 
Excursion 2, at the time of a false transient occurrence, the position of false 


transient occurrence (X,,Y,) was set at: 


X,=L-X, 
Y,=L-Y, 


Using this method, the searcher was directed away from the target in 
most cases. The results of this excursion are graphically displayed in Figures 
5.16 and 5.17. The figures clearly indicate the increased effect that false 


transient sounds have on the search effort. In particular, it is interesting to note 
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in Figure 5.17 that the curve associated with both target and false transient 
sounds has completely shifted past the curve associated with no transient 


sounds. 
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Figure 5.14. Excursion 1 - Estimated MTTD, MTBTT, and MTBFT 


TABLE V.7. EXCURSION 1 - ESTIMATED MTTD VALUES WITH 95% CI. 
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Figure 5.15. Excursion 1 - Probability of Non-Detection. 
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Figure 5.16. Excursion 2 - Estimated MTTD, MTBTT, and MTBFT 
Relationship. 


TABLE V.8. EXCURSION 2 - ESTIMATED MTTD VALUES WITH 95% CI. 
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Figure 5.17. Excursion 2 - Probability of Non-Detection. 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 


Although this thesis does not provide an all-inclusive answer to the question 
of whether or not transient detection opportunities can be utilized by a searching 
submarine to increase its chance of gaining steady-state detection on a target 
submarine, it does establish a framework, through the simulation SUBTRAN, 
that may be useful in investigating how transient sounds affect the steady-state 
search effort for a given search scenario. 

The modeling in SUBTRAN utilizes established methods. The target 
movement and transient occurrence models allow the user flexibility in defining 
the target’s characteristics. The searcher movement model is restrictive in the 
sense that the searcher has no movement options in the search pattern or 
transient prosecution. The false transient occurrence model is designed to 
supplement the search scenario with randomly occurring ambient sounds that 
could be classified as target transient sounds by the searcher. 

Major features of SUBTRAN that support it as a viable method for 
simulation of the transient environment search scenario presented in this thesis 
are: 

- Considerable flexibility is available to the user in selection of search scenario 
input parameters. 
- The simulation is designed for microcomputer implementation removing the 

"mainframe" requirement common among most submarine versus submarine 

simulation programs. 


- The simulation search scenario is easy to set-up using the input data file 
provided. 
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- The output provided by SUBTRAN is concise and clearly highlights 
important results. 


Improvements or additions SUBTRAN might include: 


- Modeling target counter-detection of the searcher. 


- Incorporating different types of target transients to include various duration 
lengths and FOMs. 


- Investigating different searcher movement patterns. 
- Including FOM target aspect and speed dependencies. 


- Modeling target and searcher movement to include speed changes. 
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APPENDIX A 
SUBTRAN USERS GUIDE 


A. SYSTEM REQUIREMENTS 

SUBTRAN requires both computer memory and disk space. The computer 
memory required is 126 Kbytes. The disk space, necessary for the two output 
files, varies with respect to the number of simulation replications performed. The 
disk space needed for a simulation run of 5000 is approximately 101 Kbytes, 
while the disk space needed for a simulation run of 1000 is approximately 27 


Kbytes. 


B. INPUT DATA FILE 

SUBTRAN requires input parameters to be entered via an input data file. 
The input data file is specially formatted to aid the user in proper positioning of 
the input parameters. A sample of the input data file is contained in Appendix 
Oy 

SUBTRAN was not coded to check for improper input parameters. Careful 
attention is required when filling out the input data file to ensure input 
parameters are within specified range limits. 

The input data file is broken up into 15 data fields. A description of the 
data requirements and allowed ranges follows: 

1. Data Field 1 

The input parameters contained in data field 1 are: random number 

seed (DSEED), length dimension of the search area (L), number of simulation 


replications (N), and the maximum search time (T,,,,). 
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DSEED - double precision variable range, [1,2147483646] 


L - real variable / range, (0,999) 
N - real variable / range, [1,5000] 
f beter - real variable / range, [0,9999] 


2. Data Field 2 


The input parameter contained in data field 2 is the number of points 


to be entered for the steady-state propagation loss curve (NP,,,). 
NP... - integer variable / range, [0,50] 


3. Data Field 3 
The input parameters contained in data field 3 are the steady-state 
propagation loss curve range values in increasing order (RN,,,) and the 
corresponding propagation loss values (P,,.). Due to the code structure for 
calculating the propagation loss curve, the values for any two sequential range 
inputs must be different. Also, the last propagation loss value must be the 


highest of all propagation loss values. 


RN... - real array of size NP,,, / range, [0,999] 


Pe - real array of size NP,,. / range, [0,999] 


4. Data Fields 4 and 5 
The input parameters contained in data fields 4 and 5 are of the same 
structure as the parameters contained in data fields 2 and 3 respectively. 
Instead of steady-state propagation loss curve, the input corresponds to the 
transient propagation loss curve. The variables in data fields 4 and 5 


corresponding to those in data fields 2 and 3 are: 
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NP,,, replaces NP... 
RN... replaces RN 


pss 


1 replaces P 


ptr pss 


5. Data Field 6 
The input parameters contained in data field 6 are the searcher’s speed 


(S,) and the target’s speed (T,). The speed parameters must never be set to zero. 


S. - real variable / range, (0,99] 


Sy, - real variable / range, (0,99] 


6. Data Field 7 
The input parameter contained in data field 7 is the number of points 


to be entered for the target’s course distribution (NP). 


NP - integer variable / range, [0,50] 


c 


7. Data Field 8 
The input parameters contained in data field 8 are the target’s course 
distribution values (C,) in increasing order and the corresponding cumulative 
probability (P,). The first value of the cumulative probability must be 0 and the 


last value of cumulative probability must be 1. 


C - real array of size NP. / range, [0,360] 


c 


P - real array of size P, / range, [0,1] 


c 


8. Data Fields 9 and 10 
The input parameters contained in data fields 9 and 10 are of the same 


structure as those contained in data fields 7 and 8 respectively. Instead of 
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target’s course distribution, the input corresponds to the distribution of the time 
between target course changes. The variables in data fields 9 and 10 that 


correspond to those in data field 7 and 8 are: 


NP,, replaces NP, 
ice replaces C,, range for T,, is [0,99] 


| 228 replaces’ P. 


9. Data Field 11 
The input parameters contained in data field 11 are the transient 
occurrence flag (F,,) and the transient occurrence rate parameter (RP,,). Transient 
occurrence is allowed when the transient occurrence flag is set to a value of 1. 
Transients are not allowed when the flag is set to a value of 0. The value of the 


transient occurrence rate parameter must never be 0. 


F.. - integer variable / range, 0 or 1 


RES - real variable / range, [.000001,99] 


10. Data Field 12 
The input parameters contained in data field 12 are of the same 
structure as those input parameters contained in data field 11. Instead of 
transient occurrence, the input corresponds to false transient occurrence. The 


variables in data field 12 that correspond to those in input data field 11 are: 


ee replaces F,, 


RP, replaces RP, 
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11. Data Field 13 
The input parameters contained in input data field 13 are the Lambda- 
Sigma Jump error process rate parameter (RP,,,) and the Lambda-Sigma Jump 


error process standard deviation (SD,,,). The rate parameter must never be 0. 


RP,,, - real variable / range, [.000001,99] 


SD,,; - real variable / range, [0,99] 


12. Data Field 14 
The input parameters contained in input data field 14 are the steady- 


state figure-of-merit (FOM,,) and the transient figure-of-merit (FOM,,). 


FOM,, - real variable / range, [0,999] 


FOM,. - real variable / range, [0,999] 


13. Data Field 15 
The input parameters contained in input data field 15 are the track 
width flag (F,,,.) and the track width (TW). The input track width is used in the 
simulation if the value of the track width flag is set at 1. The track width is set 
to the mean steady-state detection range if the value of the track width flag is 


set at 0. 


F - integer variable / range, 0 or 1 


tw 


TW - real variable / range, [0,999] 


C. PROGRAM EXECUTION 
SUBTRAN is executed by running the SUBTRAN.EXE file. The user is 
prompted to enter the name of the input data file. The name of the file may be 


any proper filename and extension of total length less than or equal to 12 
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characters. After the input data file name is entered, the input parameters are 
printed at the users terminal for verification. Simulation execution commences 
when the input parameters have been verified. Status of the simulation is 
printed at the user’s terminal every time a tenth of the simulation is completed. 
When all replications have been run, output statistics are printed at the user’s 


terminal and to the output data files. 


D. OUTPUT DATA FILES 

Two output data files are generated by SUBTRAN. The first data file, 
SUBTRAN.OUT, contains the input parameters and output statistics. The second 
data files TIME.OUT, contains the times to detection ordered from shortest to 
longest. Times to detection with a value of zero are not included. In addition, 
simulation replications that did not result in steady-state detection before T,,,, 


have time to detection set at T,,.. 
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APPENDIX B 


MAJOR VARIABLES USED IN SUBTRAN 


The major variables used in SUBTRAN are provided in an alphabetical 


listing. 


C, 


Je 
De 


FOM 
FOM,, 


Za aa at 


Z, 


NP 


c 


Course array of the targets course cumulative distribution curve 
in degrees. 


Course the target will be on after the next boundary reflection 
in degrees. 


Course of the searcher in degrees. 

Course of the target in degrees. 

Random number seed. 

Signal excess fluctuation level in dB. 

False transient occurrence flag. 

Transient occurrence flag. 

Track width calculation flag. 

Steady-state figure-of-merit in dB. 

Transient figure-of-merit in dB. 

Length dimension of the search area in nautical miles. 
Number of simulation replications. 

Number of replications in which steady-state detections occurred. 


Number of replications in which no steady-state detection 
occurred before T,,,,- 


Number of replications in which steady-state detection occurred 
at time zero. 


Number of points entered for the target’s course cumulative 
distribution curve. 
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NP 


pss 


NP 


ptr 


NP,, 


Number of points entered for the steady-state propagation loss 
curve. 


Number of points entered for the transient propagation loss 
curve. 


Number of points entered for the target’s time to course change 
cumulative distribution curve. 


Probability array for the target’s course cumulative distribution 
curve. 


Propagation loss array for the steady-state propagation loss 
curve in dB. 


Propagation loss array for the transient propagation loss curve 
in 


Probability array for the target’s time to course change 
cumulative distribution curve. 


Range array for transient detection ranges in nautical miles. 
Steady-state detection range in nautical miles. 


Range array for the steady-state propagation loss curve in 
nautical miles. 


Range array for the transient propagation loss curve in nautical 
miles. 


Rate parameter for the false transient occurrence distribution in 
hours’. 


Rate parameter for the Lambda-Sigma Jump process in hours”. 


Rate parameter for the transient occurrence distribution in 
hours”. 


Speed of the searcher in nautical miles per hour. 

Speed of the target in nautical miles per hour. 

Standard deviation for the Lambda-Sigma Jump process in dB. 
Time to the next signal excess fluctuation in hours. 

The maximum search time for each replication in hours. 


Time until the next target boundary reflection in hours. 
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| 


W 


ST 5 al 


mG 


Za 


las 


Time array for the target's time to course change cumulative 
distribution curve in hours. 


Time to the next searcher course change in hours. 
Time to the next target course change in hours. 

Time to the next false transient occurrence in hours. 
Time to the next transient occurrence in hours. 
Searcher’s track width in nautical miles. 

X component of the false transient occurrence position. 
X component of the searcher’s position. 

X component of the target’s position. 

Y component of the false transient occurrence position. 
Y component of the searcher’s position. 


Y component of the target’s position. 


65 


APPENDIX C 
EXAMPLE SUBTRAN INPUT DATA FILE 


CC C C C GeGRGee CGR CeCe ORCS eRe ee eee 


The following is the input file for SUBTRAN. 
All values should be right justified within 
their designated field. Any line with a '*’ 
or a ‘'C' in it must not be deleted, success- 
ful reading of the input data depends on 

it. If continued 4rrorsPoceunevnen - oe emp ue 
is read by the maingpuegran weensutta ne 
program documentation for the proper setup 
of this fite 


oO OF O70) S) O<Oag) OO) One 


CC CC € € € C CisG Gee eee eee 


ee) ey eee) (CG) 0G) GG)" @) 


KK KK KK KK KK K KK KK OK OK OX 
* Columns Input 

ne Soooee Se — new SSOSS KH 

* Se 15 DSEED Vad 
* ro— 3 Om L C1] 
* 33-45 N Cn] 
* 48-60 Tmax Ct] 
KKK KK KK KK KK KK KK KK OK 
Kkddddddddddddd*A*41LLLLILIL1Liliilk*xnnnnnnnnnnnnnx*x+ttctccecr eee 


me Ie 6 Fee Ce 


Zeoo 15010 1000 1000 
* 
ae KK KK KK KK KK KK KK OK KK OK 
* Columns Input * 
K ~@OF SSeS 00000 CUS Se > 
* S15 NPpss Cn] * 
KK KK KK KK KK KK KK KK OK OK 


kknnnnnnnnnnnnnk*k (CNPpss cannot exceed 50] 


4 

* 

oe KK KK KK KK KKK KK KK KK K 
* Columns IDopenbae * 
KO See re * 
* Saks RNpss(NPp) [Cr] * 
* 18-30 Ppss(NPp) [Cpl] * 
KOK KK KK KK KK KK KK OK KK OK OX 


AKOrrrrrrrrrrrr**ppppppppppppp** 
O 46 


6 70 
14 91 
18 92 
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4. KK K KK KK KK KK KK KK KK OK 
* Columns Input * 
ee SSS. lll SE * 
* SSIS NPptr Cn} x 
KKK KK KK KK KK KK KK KK K 

Peumannnnnnnnnnnk*k (€CNPptr cannot exceed 50] 

2 

x 

2. KOK KK KK KK KK KK KK KK KK 
* Columns gienbine * 
i (sane 0UC<i«‘i * 
* Sats. RNptr«(NPp) Crj x 
* ie G19, Ieie  NPp) Epil * 
KOK OK KK KK KK K KK KK KK KK XK 

AArrrrrrrrrrrrrxk*xpppppppppppppxrk* 

0 a0 
D0 120 

 S 

So. KOK OK OK KOK KK KK KK KK KK K K OK 
* Columns |ouenbie * 
aaa |060ClU, CUO * 
> SEALS) Se Ct} x 
* Om 5.0 SS Cs] * 
KOK OK KK KK KK KK KK KK KK OK OK 

mma Chee G tlt RAS SSSSSSSSSSSS KE 

2 10 

» 

7. KK K KK KK OK KK KK K KK KK K OF 
* Columns Input x 
eee Le | k 
* Sal NPc Cn] x 
KK K KK KK KK KK KK KK KK OK XK 

AKnnnnnnnnnnnnnkk (CNPc cannot exceed 50) 

6 

K 

‘om KK KKK KK KK KK KK KK KK OK 
* Columns Input * 
Tsao 8} OS * 
* Sma De: Ce (Nee? ie 
* Hen 0 fe ON Pey Cp) x 
KK KK KK KK KK KK KK KK KK 


KkKCCCCCCCCCCCCCKKppppppppppppp** 
0 


ZO 
90 
180 
280 
360 


FPF oOm Uno 
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9. KK KK KK K KK KK KK K OK XK 
* Columns Input 
* eee eee mm ee 
* Crile, Nec 


KK K Ree ee ok rere 
*kKnnnnnnnnnnnnnkk CNPtc cannot exceed 


= 
» 
10. KK K K KK K KK K KK K K OK OK 
5 Columns Input 
> 
* S= 15 swe cNint <) 
 § fea ow Pier te) 
KK K KK KK K KK KK KK K K 
AAECCCCCCCCttttkkppppppppppppprk* 
O O 
-42 i 
. 89 4 
1.43 . 6 
2.04 1 
a . 
Jee: KK KK KK KK KK KK K K OK XK 
x Columns Input 
* eeerer ere ee 8808880808000 Se 
* Soo 1h G 
k Mars), 8 Re 1g 


Ok ORD Se eK. eee. OR ee ee 
XKANNNAANNANNNAOADAAGEEYLTrrrrrrr,) 1s 


it oa 
* 
Ze KOK K KK KK KK KK K KK K XK 
* Columns Ieyerehe 
KO nnn nr ————C PS 
x Seats) eeta og 
*x 136230 UlGt ig 


KOK KK ee FO KR OR eee ee ee 
KkEnnnnnnnnnnnnnksArrrrrvrrrrrrrvrk*® 


0 6 
* 
TS. ¥ kK K KK KK K KK KK KK K 
¥. Columns Pps 
> a ee 
Ss 5-15 RPisj 
* ese. (8, SDlsj 


K KK KR K KK K KH KKK Kok x 
KArOOrrrrrrrrrrrk*ksssssssssssssxk k 
e 6 
* 
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x OK 


en 
K OK 
50} 


Ct] 


Cp] 
x X 


(m4 
fi] 
, ae, < 


Cn] 
Ges 
* * 


Toga 
Cs] 
KX 


a 


mH HM HH 


ene Fas Fat aD ap Sas 


7 toe Fe OC 


a Fe Me Mu 


KK KK KK KK KK KK K K K 
Columns ri put 


1 


e800. FOMtr 
Kaeo AK KK K KKK XK 
KKXSSSSSSSSSSSSSAACTCCECECCC CC tctctctxx 
56 90 


Je 3G 136 Je ee OG 
t 
~ 
Ol 
LB 5, 
© 
<< 
0) 
Uj 


* 
ia KK KK KKK KK KK KK K XK 


Columns Input 


* 

x 

x Se ee er 
x 3 =i Ftw 

* 18-30 TW 

x 

id 


KK KK KKK KK KK KK K OK 
fet te fe 


wre Et tf KXKSSSSSSSSSSSSSK* 


if 
0 10 
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K OF 


{s] 
[Ct] 
x 


*K xX 


{ fj 
{s] 
» Aaae 


= MH Me 


HM HM HM HM H 


APPENDIX D 
EXAMPLE SUBTRAN OUTPUT STATISTICS 


KK KK KK KK K SK 
* SUBTRAN INPUT * 
KK KK KK KK OK OK 


SIMULATION PARAMETERS 


pee 
© 
© 
© 


Replications = 
Random Seed 


| 
N 
Co 
O] 
UO) 


TRANSIENT PARAMETERS 


(1 — Yes 7 C2267 


Real Transients: 


Allowed = 1 
Rate Par. = 0.1 
False Transients: 
Allowed = O 
Rate Par. = 6.0 


SEARCH PARAMETERS 


Search Area = 10000 
Search Time = 1000 
Track Width = calc. 


Rate Par. 


steady-state Figure otf Merit) 56.00 
Range Foi nts. 0.00 Prop. Loss : 46.00 
Range Point : 6.00 Prop. Loss : 70.00 
Range Point : 14.00 Prop. Loss : 9g 
Ranges) Oint.: S00 Prop. Loss : Yange 


TRANSIENT FIGURE OF MERIT /7 PROPAGATION LOSS PARAMETERS 


Transient Pigunpesot pNer ite oo ae 
Range Point : 0.00 Prop. Loss : S008 
Range Point : 50.00 Prop. Loss : i220 7Ge 
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Course: 
Course: 
Course: 
Course: 
Course: 
Course: 


fa. OF 
ne. OF 


20; 
90. 
PSO". 
2505 
SOC: 


REE: 
Reeo. 


TARGET COURSE PARAMETERS 


Prop. 
Prob: 
Prob: 
Pro: 


Eaeab 
Prob 


a Jee 


° 
e 


« 
® 


Bc Ue 
S 
Be ioe 


% a HK 
* 0 * 


Time 
Time 
Time 
Time 
Time 


jot 
te) 
to 
Lo 
ce 


Chane? 
Change. 
Change: 
Change: 
Change: 


WITH DET BECTON TIME OF ZEKO 
WITH DETECTION TIME GREATER THAN Tmax 
CALCULATED TRACK WIDTH VALUE 


Pete teeOrtTECTION SIATISTICS 


Std. Dev. 


1/Lambda 


R-Squared 


644.19 +/-31.13 (95% 


598.44 
60255 
0.99 

Erob. of 

0.6 

Ouk7 

ones 

0.9 

eee 


gl 


« * @© @© @® © © @© @ @© @® @® @® *#® # © 8 @ © &® @ © ® @ 


O-.00,  Froeb: 0.0 
Orie aE rebp:. 0.2 
0-69 Prob: ~0. 4 
P43) Prop sae. © 
2.025 "Rroab, 1.0 

io. een Serene eee Re 0 

6) 
5.00 

(Gilley: 

Det. Byes me 

621 jal 

LOZ 

MOB 4A. 

VAIO. O 

Sali. 6 


APPENDIX E 
STEADY-STATE DETECTION DERIVATION 


As a result of the discrete event timing used in SUBTRAN, steady-state 
detection is assessed by computing the time it would take the range between the 
target and the searcher to decrease to the value of the steady-state detection 
range. The calculations used in the _ steady-state detection assessment 
computation are outlined below. 


Define, 


X, - starting X position of the target. 
Y, - starting Y position of the target. 
X(t) - X position of the target at time t. 


Y(t) - Y position of the target at time t. 


S.. .- X% component of target speed. 
Si > Y component of target speed. 
R, - Steady-state detection range. 


Xoo YeorAn(t), Y,(t),5,.,5,, - Be equivalent variables for the searcher. 


At any time t, the X and Y positions of the target and searcher can be written 


as follows; 


X(t) = X,, + tS, and Y(t) = Y,, + tS, (E.1) 


X(t) = X,, + tS,, and Y,(t) = Y., + tS (E.2) 


ey 
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RUC ne(Uem et) 1Y(t) = Y<t)y) (E.3) 
Combining equations (1), (2), (3) and setting R(t) = R,; we solve for t. 


t= B+ (B? - 4AC)°** 
ee ee (E.4) 


where, 


Poe 20.5.4 97, + S025, 5, 
Peer cee 2X, S. - 2X6, - 2X,,5,.) + [2Y.S., + 2Y,S, - 2Y..S, - 2¥.S,] 
C = X?,, + X%, - 2X.X, + Y2, + Y’, - 2Y.Y, - Re? 


The time to steady-state detection is the minimum of the positive values of 
t obtained from equation (E.4). 

Four special solution cases exist in which the time to steady-state detection 
must be investigated separately. They are: when both values of t are negative, 
when one t value is positive and the other is negative, when (B’ - 4AC) < 0, and 
when A=0. 

In the case where both t values are negative, the steady-state detection 
opportunity has already passed and detection does not occur. 

When one t value is positive and the other is negative, the range between 
the searcher’s and target’s initial positions is less than the steady-state detection 
range. As a result, detection occurs and the time to detection is 0. 

In the two cases (B’ - 4AC) < 0 and A = 0 the range between the searcher 
and the target never decreases to the steady-state detection range. As a result, 


detection does not occur. 
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APPENDIX F 


FORTRAN CODE FOR SUBTRAN 


The fortran code for SUBTRAN is divided into the following parts: 


Main Program (SUBTRAN) 

Input Subroutine (IN) 

Propagation Loss Subroutine (PLOSS) 

Target Course Subroutine (TGCR) 

Target Transient Subroutine (TGTR) 

Searcher Course Subroutine (SRCR) 

Transient Course Subroutine (TRCR) 

Steady-State Detection Subroutine (SSDET) 

Transient Detection Subroutine (TRDET) 

Lambda-Sigma Jump Subroutine (LSJ) 

Position Subroutine (POSIT) 

Boundary Reflection Subroutine (REFLCT) 

Track Subroutine (TRACK) 

Output Subroutine (OUT) 

Hold Screen Subroutine (HSCRN) 

Uniform Random Number Generator Subroutine (LRNDPC) 
Normal Random Number Generator Subroutine (LNORPC) 


Gamma Random Number Generator Subroutine (LGAMPC) 
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OO OU OO Ont Oraa® 0) O 


Ore @) OO) ©) 


es 


Ce 


® 6 @ @# @ 


* ¢ @ e e@ 


* ¢*¢ @ @ @ 


* @# @ @ «@ 


Wmem@rCur@ eG sO GNCwGrsC Ce Ce CCC CCC CEC CCC CG 

PROGRAM SUBTRAN | 
wee € Ge Cee Ge © Ce CC C.C Gee C Cee Gat CC : 
PROGRAM SUBTRAN 


x SUBTRAN requires the subroutines: FALTR, HSCRN, IN, 
Po eum Ooo me POotr  REPLCT, SRCR, SSDEI, ITGCR, 
Gg, IRACK, ERCR, IRDET, LRNDPC, ENORPC, LGAMPC. 


Declare variables and dimension arrays. 


CHARACTERKX1 ANS 

DOUBLE PRECISION DSEED 

MiP Ghe aN NEpcec eNr ptr, NPc, NPtec, Ftr, Fitr, Ftw, 

Nc, Nd, Nz, Nm, Fstc, Fsc, Fssd, Ftd, HFsc 

Rene fee Max, KRNpss (50), Pposs(S0), RNptr(50), 
pt GeO emer S CCK SO) PE 50) aeaTtc(d0>, 
Pec SO emer. Reftr, Reis}, SDIisjy FOMss, FOMtr, 
pet ceq emt 9S. YS. Xt, Yt, UC2), Ge, ITCt, TIt, 
Ce eles. At, Yf, Tit, Bist. LbicH, Rssde Cr, Ir, 
jimonves losa. la<5000);, HxXs, HYs, HCs, HICs, Hih 


Read input data from the user-defined input data file. 


CALL IN(CDSEED, L, N, Tmax, NPpss, RNpss, Ppss, NPptr, 
Mime Dinas t, OS, NrC,me@c, Pc, NPtc, ttc, 
ee meron het kftr> RePftr, RPls}, SDisj, 
FOMss, FOMtr, Ftw, TW) 


roe INPUL Parameters are COrrect. 


WRITE CX, 14) 

READ (X, 15) ANS 

ME CANS ,NB. °’Y’) .ANDoe (ANS .NE. * y' )O THEN 
WRITE CX, 24) 
CALL HSCRN 
SOTeG 7 oO 

ENDL 

WRITE CX, 34) 


If the track width is not specified by the user 
(Ftw=0) then compute the searcher’s mean track width 
(TW) in nautical miles using the steady state figure 
of merit (FOMss) and the propagation loss curve. 


Pesch tw ee 0) THEN 
CALL PLOSS(NPpss, RNpss, Ppss, FOMss, Rd) 
TW=2%*Rd (1) 

PND Er 


15 


CO) 0) © @aerO 


RPOOADNDDAAD 


CG a@ 1). G.-C) 


COs Geel) 


COG) eG) G) 


ee © @¢ @ @ 


eo @ @ @ 


oe @ @ @ 


« @¢ @ @ @ 


©. ey 8.6 Ce 


Start counter for number of simulation replicatiGne 
(Nc). Initialize detection counter (Nd), detection at 
time zero counter (Nz), and detection not before Tmax 
counter (Nm. 


Nc=1 
Nda=0 
Nz=0 
Nm=0O 


Assign the initial value of the simulation time (T) in 
hours to zero. This is the starting place for 
successive runs. Set the searcher'’s transient course 
flag (Fstc) to zero. 


T=0 
Fstc=0 


set the initial position of the searcher (Xs,Ys) in 
the lower left-hand corner of the search area one-half 
a track width’s distance from the edges of the search 
area. 


XS=TW/2 
Ys=TW/2 


set the dmitial position oi the targee xt, a 
uniformly distributed in the search area. 


CALE LRNDPC(DSEED IU.) 2 
Xt=LAU C1) 
Yt=LxAU Ce) 


set initial target parameters of course in degrees 
(Ct), time to next course change in hours (TCt), time 
to next transient in hours (TTt). Target’s speed is 
assumed constant. If the target’s transient flag (Ftr) 
is i, transients are allowed in the simulation. 


CALL TGCKCDSEED, NPc, Cc, Pc, NPtc, VTtc, Pte nie cee 
[FF CPt eer ODEN 
CALL TGITRCDSEEBD Bert tn.) Le 
ELSE 
TTt=2kTmax 
END IF 


set the searcher’s course flag (Fsc) to 1. 

Set initial searcher parameters of course in degrees 
(Cs), time to next course change in hours (TCs). 
Searcher’s speed is assumed constant. 
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s&s ¢ © # @ 


eG) Oro 


ces € «6 


es ¢ © # ¢ 


¢ @¢ @® @ 


CO) (Qi) (GVO MO IS Oe {Geile ie. 


i 


Fsc=1 
CALL SRCR(Fsc, Xs, Ys, Ss, TW, L, Cs, TCs) 


If false transient detections are allowed (Fftr=1), 
Pyemmoet the tmitial false transient parameters of 
bearing to the transient in degrees (Cf), time to next 
false transient (TTT). 


ie Cettr . bOs 1>THEN 

Clie ratiReepDSBPED, RPftr, L, Xf, Yi, TTt) 
Bie Se 

TT£f=2*kTmax 
END IF 


set initial signal error level (E1sj)») and time time to 
next error level fluctuation (Tlisj). 


Web seeborrey wert 1S), SPs}, Elsj;, Tisj> 
Sempre ce Steaa,y Stave derection range. 


FOM=FOMSStE1sj 
CALL PLOSS(NPpss, RNpss, Ppss, FOM, Rd) 
Rssd=Rd (1) 


Compute the time in hours until the target must be 
reflected off the search area boundry (Tr) if present 
eotinsewencdeSpecaduane maintained. Also compute the 
ecumcemolvere er ection (Cr). 


wire wen cd. Cl, St, At, Yt, Cr, Tr) 
Computes the winimum time <Imin) until the next event. 
mito Mi Ccimax ioe let, F1t, TCs, T1f,1Tlsj, Tr) 


Check to see if steady state detection occurs before 
iit So Stop simulation run (Fssd=1). Time to 
steady state detection (Tssd) is computed if detection 
occurs. 


@abibsSDET CXL, »¥t, ot, Ct, Xs, Ys, Ss, Cs, Rssd, Fssd, 
Tssd) 
DPew@crecume@ 0) . AND. <Tssa. LE. Tmin) > THEN 
be Clee ee Oo, AND] (Issd .EQ. 0)> THEN 
Nz=Nz+1 
GETO 50 
eee SE 
GOTO 40 
END IF 
END IF 


ve 


GIOee OG OO) 0) 0) 1) 0 O © 


O.@)@ 


COG): 


* @¢ © © @ 


6 © #® @ 6 


Calculate the new searcher and target positions based 
On” lms, 


CALL POSIT(CImin, At, Yt, €t, St. 4s Scere oe 


Decrement the time values and check to see what the 
next event is 


P=itimin 

per= [er imin 
tit—-1 te i min 
TCs=TCs-Tmin 
1T£=TEf=ilmin 
Tisj=Tisj-Tmin 
ay = fee mim 


Route to various updates based on which of the 
following events is due to occur. 


1. Maximum search time is up. [UT >= Sims 
ce. Target course change. [TCt = 0] 
3. Target transient. (TTt = 0] 
4. Searcher course change. (TCs = OQ) 
5. False transient. CTT£ = O) 
6, Shenal error fluctuation. [(Tlisj = QO] 
7. Target boundry reflection. (Tr = 0] 
Update 1. 
if CT GE) imax) ieee 
Nm=Nm+1 
GOTO 50 
EB Dae 
Update éd. 


1p CTCEe ee... Oo pal 
CALL TGCRCDSEED, NPc, Cc, Pc, NEte, icc we lc eee 
Tore 
COLOR 
END IF 


Update 3. 


Le Ghee 2G) Oo THER 

GREE TGTR (DSEHED eRe tr. 1 te 
FOM=FOMtrt+Elsj 
CALL PLOSS (NPptr, RNptr, (ptr eerreca 
CALL TRDET (Xt; Yt, Xs, Yo. Rape ce 
IF(Ftd .NE. 1)THEN 

GOTO 30 
ELSE 

Fstc=1 
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CALL 


coTo 
END IF 
END IF 


Update 4. 


IF (TCs 
1PCFste 
CALL 
GOTO 

ELSE 
CALL 


GOT 
END IF 
END IF 
Update 5. 


ee CT ed. 


. EQ. 


. EQ. 


TRER (hse, 
HFsc, 
TCs) 


FPSstc, 
HXs, 


o Ger 
HYs, 


Vtneeeo. YS, 
HCs, H1ICs, 


Sen L, TW, 
Hin, es; 


30 


0) THEN 

pe QeeO) THEN 
SRER (PSC, Xs, 
SO, 


Ys, Ss, TW, L, Cs, TCs) 


DPRER CPSC, 
HFsc, 
TCs) 


Fstc, 
HxXs, 


boas 
HYs, 


Yt AS, 5S, 
HCs, HTCs, 


Ss, L, TVW, 
Hen ee, 


30 


O) THEN 


FOM=FOMtr+Elsj 


CALL PLOSS(NPptr, 
CALL TRDET(Xf, 
IF(Ftd . 


CALL 
GOTO 
ELSE 


Fstc= 


CALL 


CALL 
GOTO 
END IF 
END LF 


Update 6. 


IF(T1sj 


. EQ. 
CALL LSJ (DSEED, RPlsj, SDlisj, 


RNP Cr weep cr, 
Xs, Ys, Rd, 


FOM, Rd) 
ere lah aol) 
NE oA N 

athe CDom ED. Reitr, La, 


30 


Ro et nee ele) 


1 

Re (ese. 
HFsc, 
TCs) 

PA GDSEED. RPitr, L., 

30 


bSce. 
Hxs, 


et. 
HYs, 


et. 
HCs, 


Xs, YS, 
HTCs, 


Ss. vey TW. 
HTh, Cs, 
Xf, 


Meier) 


O) THEN 


Bikoq, lisj) 


FOM=FOMsSstElsj 


CALL PLOSS(NPpss, 


RNpss, Ppss, FOM, Rd) 


Rssd=Rada (1) 


CeTe 30 
PND Vir 


Update 7. 


a 


LER re bo. Oo tae 


Cr-er 
COTB2 26 
END IF 
C 
Ores: Store times to detection (Td). Update detection 
C counter (Nd). Update replication counter (Ne). Keuge 
C for next replication run if required. Print te) gee 
C screen a message every 100 replications. 
C 
40 Nd=Nd+1 
Td (Nd) =T+Tssd 
5/0) Ne = Neen 
IF CCINT< (Ne-1)7 (N/10))04(N710)> . EQ. (Nea teen 
WRITE CK, 44) (Nc-1) 
SIRBB IS 
[FPF (Ne TES Sh COrTeaie 
C 
Cie Compute the output statistics and print Ver the ioe 
C terminal and output file. 
C 
CALL OUTC(Nd, Nz, Nm, Ftw, TW, Tmax, Td) 
C 
Cr eee Print a closing message to the user’s terminal. 
C 
WRITE CA, 54) 
C 
Cae Read format statements. 
C 
IS) FORMAT (A) 
C 
Cnn nae Write format statements. 
C 
14 FORMAT (23 (/),1X,’ Are the input parameters correct?’, 
I. OCR Eer ge 
24 FORMAT <(23(797> 27%, "EXECUTION TERMINATE DSEY (USER 
il Fane 3 GMP) 5) IG) 
34 BORMA LCZOK7 Dele Ee Cue alg ee ae 
44 FORMAT (5X,14,' replications complete.') 
04 RORMATCZe07 ), Sako G tee Zook SUBTRAN x’ / 
i yO Kae ey a7 
il 30X,’ Simulation Complete’ ,9¢/)) 
60 END 
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M~mememescnee © CC CC Ce@me ec CC Cc cc cccCcccmcC cc 
C C 
C SUBROUTINE IN C 
C | e 
P~ememene ween CCC SC Ceeme Cc C Cree ccccecccec 
C 

SUBROUTINE INCDSEED, L, WN, Tmax, NPpss, RNpss, Ppss, 

1 Veet erN@Gwt Pile, oF, oS, NPe, Cc, Pe, 

il Neve | ten we, Fire krtr, Frtr, RP£tr, 

3 RPlsj, SD1isj, FOMss, FOMtr, Ftw, TW) 
C 
C Variables passed to INPUT are: none 
Cc Variables returned by INPUT are: all in Subroutines 
C Cart List 
Cc 
C ASIN requires the subroutine HSCRN. 
Cc 
oo ae Define variables. Dimension arrays. 
C 

CHARACTER RAI2 Peet N, FILEOUT 

DOUBLE PRECISION DSEED 

Vio we Pos menh pum, NPC, NPtc, tr, Ritr, Ftw, 

1 ITHOLD 

Rei, May DSS (0), PpsstS0), RNptr dod), 

1 Piper (oO owmmros, Ce (D045 Pe (504° Tte (50), 

ih MUG S07 iaem ert. RPES lye oDls}], FOMss, FOMtr, 

1 TW 
C 
oo Write to the user’s terminal an intoruction message. 
C 

WRITE CA, 14) 

CALL HSCRN 
C 
Co oe Read the name of the input data file from the user's 
terminal. Assign 
C the name of the output file to be used to store the 
moout data. 
C 

WRITE CK, 24) 

Peni oor TLE] N 

Pe BOUT="SUBTRAN. OUT’ 
C 
Ca Open the user-defined input data file and the assigned 
C GButput data file. 
Cc 

Siveeero, PILE=hiLe iN) 

SCEPENCZO, Er ILE-FILEOUT) 
C 
oO Read the 15 data fields from the input data file. 
C 
Meee Yate field 1. —DSEED,L,N, Tmax] 
Cc 
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O 


© 


Orr @ 


Cio, Geno © 


READ (10,25) PeBEe A, Nt mex 
Datawttelad 2. UNE pss 
READ (10, 35) NPpss 
Data field 3. (RNpss(NPp), Ppss(NPp)] 
DO 10 [=1,NPpss 

READ (10, 45) RNpss(I), Ppss(I) 
CONTINUE 
Data field 4. (€NPptr} 
READ (EOj SS) Ni ere 
Data field St Nprer (ee per (NEE 
DO'2Z20 [=1, NEper 

READ CIO, 45) RN pti el), hepa 
CONTINUE 
Data field , 7 TSt, cs) 
READ CLO. D2)s5t,05=> 
Data field * UNE! 
READ ClO, co Nee 
Data field Ga. f Ge Nee)” Pe (Nc)? 
DO 30 [=1,NPc 

READ ClO. 49 Ge Cl are Cy) 
CONTINUE 
Data field 9. [(NPtc] 
READ CTO; oo NEte 
Data field 10. Vitte (NEteoy rte tir ve a 
DO 40 [=1,NPtc 

READ CIO 45071 teq 7 ve op 
CONTINUE 


Data ftelad wie Baw kr tre) 


READ (10, 65) Ptr skiers 
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oe) 


OC) 2 Ora 


Go) 


Er) OVO) 0). Ore sOn@ OG) 


* ee © 8 ¢ 


Warcaeceela i2. ( Fitr, RPfitr] 

Real 0, 65) Fi tw KPLtr 

Aad Storer tSj, ols] i 

READ (IO, O30 RFis}, SDisj 

Data field 14. (€FOMss, FONMtr] 

READ (10,55) FOMss, FOMtr 

Data field 15. [Ftw, TW] 

READ CLO;O5> Fr tw, TW 

ELineumeoeunenusSen S Cerminal and the Gutput data file 
a verification of the input. There are three screens 


of information. 
1. Miscellaneous data. 


Zee npOopagation loss Curve —- steady state 
(points). 

Seaapacatton GSS Curve — transient 
Coes GS . 


4. Target’s course and time to course 
change curves (points). 


Screen 1. 


WRITECX, 4) 
WRITE CX, 34) 
WRITE C20, 34) 
WRITEC*; 44)N, (L**2), DSEED, Tmax 
WRITE (20, 440N, (L*¥*2), DSEED, Tmax 
IFCFtw . EQ. 1)THEN 
WRITE (CK, 54) TW 
WRITE (20,54) TW 
ELSE 
WRITE CX, 64) 
WRITE (20, 64) 
END IF 
preteew 4) Relics, SDIsi, Ptr, RPtr,Fitr, Ss, RPftr,st 
Wie C20ma4o RPiSiacDIS|,Ftr, kKPtr, Fitr, Ss, RPftr,st 
CALL HSCRN 


Screen @2. 


WRITE(*, 4) 
WRITE(*, 84) FOMss 
WRITE(20, 84) FOMss 
IHOLD=0 

DO 50 I=1,NPpss 
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50 


60 


70 


80 


GG) ©) 


LHOLD= Piotr 
WRITE CK, 94) RNpss(I), Ppss (I) 
WRITE (20, 94) RNpss(I), Ppss(¢I) 
IPCC. EQ3 815) “OR Clg SOO ee 
THEN 

WRITE CK, kX) 

WRITE CX, *) 

CALL HSCRN 

WRITECK, 104) 

ITHOLD=0 
END IF 

CONTINUE 


IFCC]l . NE. 16)..0R. (1 (NES si) 20. Fees 


DG GOPl=1, 1921 8SrEP 
WRITE CK, *) 
CONTINUE 
CALL HSCRN 
JNO) Iie 
WRITEC2Z07 154) 


screen 3. 


WRITE CX, 4) 
WRITE (xX, 86) FOMtr 
WRITE(20, 86) FOMtr 
IHOLD=0 | 
LB) YAO i asl Nshone 1 
IHOLD=1HOLD+1 
WRITE (*, 94) RNptr (1), Pptr (1) 
WRITE(20, 94) RNptr (1), Pptr (1) 
[FCC] .EQ. 15) 2 OR, (Cl sEQ. 30) aera! 
THEN 
WRITE (x, *) 
WRITE (X, *) 
CALL HSCRN 
WRITE (X*, 104) 
IHOLD=0 
END IF 
CONTINUE 


LF<< 1 2NE.. 16> .0Ok. «1 .NE. 81) Vek. (lee 


DO 80 I1=1,19-IHOLD 
WRITE CX, k) 
CONTINUE 
CALL HSCRN 
END IF 
WretiE<20) 154) 


Screen 4. 


WRITE C#, 4) 
WRITE CA. 1 14> 
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-HQ. 4599 


46)) THEN 


.EQ. 450%) 


46)) THEN 


100 


90 


oO 


24 


PRP rr 


wnt ee O14) 

JHOLD=0 

DO 90 I=1, MAX (NPc, NPtc) 
IHOLD=IHOLD+1 


ie GCGl eaereeNre) . AND. <1 


Ge NRtc) ) THEN 


whe Gt 2c Cl) ere CI) ite tl), Pte (I> 
Wimme (Zen tee) Cet lore, ttc Cl), Pte I> 


GOTO 100 
JORDI Is 
PE ci. GlowNP econ EN 


Wil ib Gel S49 ee Ch), Pe tl) 
WRITE C207, 345 Cetl), Pc< 1) 


Cero 100 
END 1 
Decl sei NPeo THEN 


WreltrE<x*, 1440 TEeC1), Ptc cl) 
Whe te Geo 144 Tte cl), Pte cl) 


GOTO 100 

PND i 

PG. -e aster OR. Cl 
WRITE CK, *) 
CALL HSCRN 
WRITE (X*, 104) 
IHOLD=0 

BND ik 

CONTINUE 


ECC Ge NE. 19) 2OR, Cl «NE. 


Dey lO eee] 1 HOLD 
WREde Cx 
CONTINUE 
CALL HSCRN 
END IF 
WRITE (ZO, 1954) 


Read format statements. 


FORMAT (A) 


-EQ. 36)) THEN 


Sf THEN 


mre Gz o.67 2 ls 0.2%, P1300, 2x, 113,2xX, F13.90) 


FORMAT (7 (/), 2X, 113,8¢/)) 
FORMAT (2X, F13. 0, 2X, F13. 0) 


mene (O7), 2X, PLS. 0,2X, F13.0) 


mekMAdeCSa (7), 2%, 113,24, Fis.0) 


Write format statements. 


FORMAT (23 (7) ) 


FORMAT (23 (7) ,20X,’Submarine Transient Detection ’ 
Comte cen O4/ > Son, OC KX ?)/SeLK, 
ae SUBTRAN Nea oeree Gt), O67), 29K, 
Memeo aaGmiaabe tied .~USN' 73ox,'26 July ’, 


moe O67 >> 


Hm 2s7 1X, EPNIER THe NAME OF THE INPUT DATA 
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9 


i! an Te Ee 


34 FORMAT (29%, 100’ * ')/29K,' x  SUBIRAN SINE Tee 
1 29K; 100s DD 
44 FORMAT (//T9, 'SIMULATION PARAMETERS’ ,T50,’SEARCH ’°, 
1 *PARAMETERS'/ 19,21 ©@ = 03150 eG 
1 T9,'’Replications = *,16,150,’ Search Area —aae 
Jt FO.2/T9, ' Random Seed = ',F9.2,T50O, 
it ‘Search Time = ',F9.2) 
54 FORMAT (T50,' Track Width = ',F9.2) 
64 FORMAT <(T50,’ Track Width —=sealte.) 
Fee. FORMAT (TS, ’ TRANSIENT PARAMETERS’ /T9,20¢’'-'), T50O, 
1 "SIGNAL ERROR ” .” PARAMETERS’ 7 iter 
i ’(1 = Yes 7 0 — Ned’ 7150,286 — 77 1c 
Hi "Rate Par. = ',F5.2/T9,’Real Transients: ',T50O, 
I ‘Std. Dev. = ',F5.2/T10,'Allowed = ', |27 ier 
1 *Rate Par. = ',F7.4,T50,'SPEED PARAMETERS’ /TO9, 
1 ‘False Transients: ’ ,1T50,16C¢’-')/T10, 
1 ’Allowed = *',12,1T50,’S. Speed = *,FS.27 iis 
1 "Rate Par. = ',F?7.4,1T50,’'’T. Speed = ’, FS72% 
1 4-7) ) 
84 FORMAT (12X, ’STEADY-STATE FIGURE OF MERIT / ', 
it ’PROPAGATION LOSS PARAMETERS’ /12X, 
it 99('=')//2e22ek, 'Steady-state FPigureseor ae, 
it "Merit = *~. FO.27) 
86 FORMAT (13X, ' TRANSIENT FIGURE OF MERIT / PROPAGATION '°, 
it *LOSS *,’ PARAMETERS’ 713X%,50('—-°)77 “2a 
it ‘Transient Figure of *,’Merit = *», FG7270 
O4 FORMAT (fis, Kange Point <2 3) koO72 oo, 
it DPSOpawelesSs -)'), hOx2z) 


104 PORTA TG (7 y) 
114 FORMAT (28X,' TARGET COURSE PARAMETERS' /28X, 24 ¢'—')/7/7) 


124 FORMAT (5X, ’Course: ’,F4.0,3X%,'° Prob: ©) h4-2. 42. 
1 ‘Time to Change: *,F6.2,3%,° Pre@bige seer 
134 FORMAT (5X, ’Course:? ° , F4. 0, 3X," Frab:. | F442) 
144 FORMAT (T45,’*’ Time to Change: °,F6.2,3xX,"Prab. | fie 
E54 FORMAT (3 (7) ) 
END 
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CO mr) Or) ant) C) 


20 
30 


Cem ememeuoeCweCuC ECs C-GeG CCC CCC CCC CCC 


SUBROUTINE PLOSS 


Bo Oe 


PIC eCNeOseEOnClGUON@rOnCuc CuG CC CCC C CCC 
SUBROUTINE PLOSS(NPp, Rp, Pp, FOM, Rd) 


Variables passed to PLOSS are: NPp, Rp(50), Pp(«50), 
FOM 
Variables returned by PLOSS are: Rd(7) 


Declare variables. Dimension arrays. 


INTEGER NPp, COUNT 
REAL Rp(50), Pp(50), FOM, Rd(7), SLOPE, Rint 


Initially assign the detection ranges, Rd, to zero. 


DOM UO =i, 7 
Rd (1)=0 
CONTINUE 


Compute the intersection of the propagation loss curve 
and the figure of merit. 


COUNT=1 
DO 20 I=1,NPp-1 
SLOPE=(Pp(I+1)-Pp(1))/(Rp(It+1)-Rp (1) 
Rint=((FOM-Pp(1))/SLOPE)+Rp (1) 
IF((Rint .GE. Rp(I)) .AND. (Rint .LE. Rp(I+1)))THEN 
Rd (COUNT) =Rint 
COUNT=COUNT+1 
IF(COUNT .EQ. 8)GOTO 30 
END IF 
CONTINUE 
END 
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CC 
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1 


*s t§ @¢ @ «6 


ON) OOO Ore. © 


ee 8 © @ @ 


OrQ@i@e) ©) 


eo e -¢ ‘e 8 


GO). ©) 


30 
40 


ccccccceccecec Cec € € CC € € © Cle eee 
C 

SUBROUTINE TGCR Cc 

| C 

CCCCcCCCcCCcecC Ce CC CCC Ge CCC Cee 


SUBROUTINE TGCR(CDSEED, NPc, Cc, Pe, NPte, Tte,) bees 
Croce 


Variables passed to TGCR are: DSEED, NPc, Cc, Pc, 
NPte,. 1 te, eve 
Variables returned by TGCR are: Ct, TCt 


x TGCR requires the subroutine LNORPC. 
Declare variables. Dimension arrays. 


DOUBLE PRECISION DSEED 

INTEGER WNP ceive 

REAL Cc(50), PcetS0), Ttct30), Pte(S0), Ct, [Ci tee 
SLOPE 


Generate two Uniform(0,1) random variables to be used 
in extracting the new course and time to course 
change. 


CALL LRNDPC(DSEED, U, 2) 
Extract the new course value (Ct). 


DO 207 )=1, Nees 
SLOPE=<(Pe CIil+i)-Pe ti) 7 GeGit eve Gee 
IFCCUC1) .GE. Pecl))>) .AND. (UCI) EE, Pe <i ieee 
Ct=( (U Cl =Fet)) 7 SUOr Ee) eee 
GOL@? Zo 
END IF 
CONTINUE 


Extract the new time to course change (TCt). 
DO{307 1=1,,NEteet 


SEGPE=(Pte<Il+1)=Pte Clo, CiteGit) —~fre eee 
ifc<«UCe) .GE. Ptetl)>) .AND. (U(2) 7 LED Pic Giaaeee 


THEN 
TCe={ CU <2) -Fte Cl)? Sreole tive ae 
GOTO 40 
END IF 
CONTINUE 
END 
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CG) 2G) 


Cee 


ae 


* #© &® 8 28 


Meememeor ee remo eC ee CC CC CCC CC CC CC 

SUBROUTINE TGTR | ‘ 
CGememere = © C Cee Ce CcCcececccceoeccececc : 
SULEROUMINE  TGER(DSEED, RPtr, TTt> 


Vartables passed to TGTR are: DSEED, RPtr 
Varlables returned by TGTR are: TIt 


x TGTR requires the subroutine LGAMPC. 
Declare variables. 


DoUSCE PRECISION DSEED 
REAL RPtr, TTt,& 


Teper e tne slime tO Next transient «1Tt). The time 
between transients is exponentially distributed with 
team -~sl/RPtr in hours. 


GAME LGAMNPC CDSEED, E,1,1.) 


We —-E/7RPeY 
END 
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CCC CCCCcCCcCeceCC CC C © CC Ce eC Cee ea 


SUBROUTINE SRCR 


aD) 12) Qe. 


CCCcCCCCCcCCcececcecCcCc Cc C C € € CG <€ Clee 
SUBROUTINE SRCR(C, Xs, Ys, Ss, TW, L, Cs, TCs) 


Variables passed to SRCR are: C, Xs, Ys, Ss, TW, L 
Variables returned by SRCR are: Cs, TCs 


x The variable C is generated by SRCR when it is first 
run and is passed back to SRCR on successive runs. 


Declare variables. 


INTEGERSSC 
REAL Xs, Ys, Ss, TW, L, Cs, TCs, Pl 


Assign the value of Pl. 
PI=3.1415926 
Set up the program routing to control the course. 


LP CCE S2EQS 1) SOR eee EG 5) CeTe To 
LE«<G, 7 BQ. 205 On we ©. Eo, 6) GOle Zo 
IFCCE SEQ: S30 OR. KG Eo. 7) COteeso 
LIFPCCC EQ. 4093 Ch (eC 7 OOo) Geese 


Searcher’s course - north. 


Cs=0 

TCs=(ei0 7455 

le (CC. EG 1) THEN 
C=2e 

Biss 
C=6 

END IF 

GOTO 50 


searcher’s course - east or west (following north). 


LER re howe THEN 

Cs=90 

1P CxS ~GEw (i= <1 SkiWD) THEN 
TCs=¢€ (CL- (TW/2) )-Xs)/Ss 
C=7 

ELSE 
TCs=TW/Ss 
C=3 

END IF 


90 


0 


ELSE 


Cs=270 
lF(xs . LE. (Cl 2eAy)) THEN 
TCs= (Xs- (TW/2))/7Ss 
C=3 
ELSE 
TCs=TW/Ss 
C=7 
END IF 
END IF 
GOTO 50 
Searcher's course - south. 
Cs=180 
TCs= (L-TW)/Ss 
Pace 1 EG: 3) THEN 
C=4 
JE psy) 
C=8 
END IF 
COt®. 50 
Searcher's course - east or west 
rece . EO. 4) THEN 
Cs=90 
PraGtes ©. Gl. (l(a * Il WwW) THEN 
TCs=(€ (L- (TW/42))-Xs)/Ss 
C=5 
BESe 
TCs=TW/Ss 
C=1 
END IF 
ELSE 
Cs=270 
Peets 2EE. Ci. S4TW)) THEN 
TCs= (Xs- (TW/2))/78Ss 
G=1 
Bioe 
TCs=TW/Ss 
C=5 
END IF 
END IF 
GOTO 50 
END 
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OO iGr iG) Ol) 


CCCCCCcCcCcCcChUCCECCC CCC CC € CC Clee eee 


SUBROUTINE rer 


Ccccccccccrecrc CCC CC Ce € CC Cie eee 


C 
C 
C 
C 
C 


SUBROUTINE TRCRCC, Cl, Kt, Yt, Xs, Ys, Ss, L, Wie 
HXs, HYs, HCs, HTCs, HTh, Cs, TCs) 


Variables passed to TRCK are: C, Cl, Xl, Yue 


5S,0 Uae 


Variables returned by TRCR are: Cs, TCs 


* 


* 


The variables HC, HXs, HYs, HCs, HTCs, HTh are 
generated by TRCR when it is first run and are 
passed back to TRCR on successive runs. 


TRCR requires the subroutine TRACK. 


Declare variables. 


INTEGER VG= el. He 
REAL Xt, Yt, Xs, Ys, Ss, L, TW, HXs, HYs, HCs, Hig 


hit, -¢s, f¢s 


.eoe.. DET up routing depending on Cl value as follows: 


Ys, 


1; trancsilempumuaciaeder cater CCi=1) 
e. Return back along the transient course. [C1=2] 
3 


Return to "on hold” steady state search 


positon: (C1=3]) 
4. Resume steady state search. (C1=4) 


hE Cel (2651) GOTO so 
FF (Cli Eel 22 GOTO s20 
DECC I FO] SU eee 30 
LE CCl” ~EO:- 476010 46 


Route 1. 


VG eee nO) en 


HXs=Xs 
HYS=Ys 
HCs=Cs 
HTCs=TCs 
HC=C 

C=9 


END IF 
CAL DPERACK (Xt. Yi. Xs. You Ss) leo. fesp 
TCs=TCs- (TW/ (2*Ss) ) 

Hih=ne= 

Cle 

GOTO! 50 
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WOOD 


ayo) 


Route 2. 


Peace oar tee 160. 7HEN 
Cs=Cst+180 

Peo 
Cs=Cs-180 

END IF 

TCs=HTh 

C1i=3 

GOTO 50 


Route 3. 
CALL TRACK (HxXs, 


C1=4 
COO. 50 


Route 4. 


C=HC 
Cs=HCs 
TCS=HTCs 
C1=0 

END 


HYs, 
TCS=SQRT (¢€ CHYS-YS) KX2)+¢€ CHXS-XS) KK2))/7Ss 


Xs, 
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CCCCcCCCCCC CC CC € CC € CC CeCe ee 
@ 

SUBROUTINE SSDET c 

c 

CCCcCCcCcCCCcceceC CCC C€C OC CC Che ere. eee 


SUBROUTINE SSDETCXt, Yt, St, Ct, Xs, Ys, Ss, Cs, Rssda, 
Ci, Tssd) 


Variables passed to SSDET are: Xt, Yt, St, Ct, Xs, Ys, 
Ss, Cs, Rssd 
Variables returned by SSDET are: Cl, Tssd 


Declare variables. 


INTEGER Gi 

REAL Xt, Yt, St, Ct, Xs, Ys, Ss, Cs, Rssd, Issa 
HOLD, Stx, Sty, Ssx, Ssy, A, B, ©, Ax, "Sx ew 
Ay, By, Cy, Tssda, Tssdb 


Assign the value of Pl. 
PI=3.1415926 


Check to see if the searcher is already within steady 
state detection range. 


HOLD=SQRT ((€ (XS-Xt) KK2)04+ CCYS-Yt) KkK2) ) 
IF CHOLD SES] Rkssao THe 

Tssdad=0 

Ci=1 

GOTO 20 
END IF 


Compute the target’s x speed component (Stx) and y 
speed component (Sty). 


Stx=StxSINcCtxPI/180) 
St y=StxCOsS (Ct *P17180) 


Compute the searcher’s x speed component (Ssx) and y 
Speed component (Ssy). 


SSx=Ss*SIN(Cs*P1I/180) 
Ssy=Ss*COS (Cs*P1/7180) 


Compute the constants to be used in the quadratic 
equation for solving the time to steady state 
detection (if it occurs). 


Ax= (Ssx*4#eC)+ (St xx*2) — (2xSsx*Stx) 
Bx= (2kXSkSSxXx) + COKXCKXStx) — (CSKAXSKSEKX) — CS™SKKXCKASSxX) 
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OOO) Sore Ose) On) CrOrOe® DIAS. (YOO ©)"C) Ore 


Se Bi 


* ¢ @© # @ 


Cx= (XSKKS) + CXC KK] — COKXSKXKL) 

Ay= (Ssy*xe) + (Sty*x2)- (C2xSsyx*Sty) 

By= (2*4YSkS5SsSy) + (C2kYtCHKSty) — C2kXYSkSty)d—- CexYtxkSsy>) 
Cy= (YS**2) + CYtCKk2)-— (CAXYSKYt) 

A=AxtAy 

B=Bx+By 

C=Cx+Cy- (Rssd*xk2) 


G@heenm tne veonmarttions that yield a solution that is not 
heal Or Linite, 

1. A=0 

2. (Bxx2)-— C4XAXKC)<O 


Caonaittion 1. 
PR CAS EO. "O0) GOTO 10 
Con@satlion 2. 


HOLD= (B**2) —- (4X AC) 
PacHOlD ~LT. 70>GOTG 10 


Conpiuce Lne values £Or the Eime Lo Steady state 
aetection <(Tssd>. 


Tssda= (-BtSQRT CHOLD))/7 C2X*A) 
Tssdb= (-B-SQRT (HOLD) )7 (C2*A)D 


Check for negative detection times. 
1. Both Tssda and Tssdb negative. 
Cee fSsaa specittive, 1Sssdb negative. 
oe issdaenegarive.es (SSdbypestiive. 
4. Both Tssda and Issdb positive. 


Case 1. 
eeackeaccdiame FF. OF .AND. “Tssdb .h1. “0))>GO!IO 10 
Case 2. 


Megs caame Ge. O) ~AND. <Tssdb .LIT. 0O)) THEN 
Tssd=Tssda 
Ci=1 

END IF 


Case 3. 


fPpmecreeda .L1, 0) 2AND. (Tssdb .GE. 0)) THEN 
Tssd=Tssdb 
C1=1 

END IF 
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me) 
Zo 


Case 4. 


IF«(Tssda .GE ROD “AND Clssdbpaace, 


Tssd=MIN (Tssda, Tssdb) 
C1=1 

END IF 

GOTO 20 

C1=0 

END 
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Pee memememenere Clerenme © 1nG C CC CC GE CC CCC 

SUBROUTINE TRDET : 
Semone ree CoCr eee CaGec © CC CCC EGC Cc CC 
SURE OUTINGSIRDET CGtjeyt, ASseers, RdywaG 


Variables passed to TRDET are: Xt, Yt, Xs, Ys, RdC7) 
Variables returned by TRDET are: C 


Declare variables. Dimension arrays. 


INTEGER C 
MEM, LC, AS, LS, RdC7), HOLD 


Compute the current range between the searcher and 
target (HOLD). 


HOLD=SQRT ( ( (XS~-Xt) *k2)+ CCYS-Yt) KK2)) 


Check to see if the current range is within transient 
detection range. The code (C) for a detection is 1, 
and the code for no detection is 0. Two cases exist 
that require checking: 

iteeoLneem path cetection: 

ey Cenvergence zone detection (3 max). 


Case 1. 


Co Dane. haecl)) THEN 
C=1 
Gato 20 

END IF 


Case 2. 


POO. l=2,0, 2 
Neeenelle. GE.) kalo) AND. “CHOLD”.LE. Rdi(Il+1)))THEN 
C=1 
SOT 2.0 
END TE 
CONTINUE 


Indicate that no detection occurred. 
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CcCCCCCCC CCC € CC CC © CG ee ere eee 
SUBROU] INEMESs3 : 

CCCcCC CC CEC CC CC Ce © CC 6 Ce een en ee : 

SUBROUTINE LSJ (DSEED, RP1lsi je sDlsi Elsie 


Variables passed to LSJ are: DSEED, RP1lsj,SDlisj 
Variables returned by LSJ are: Elsj,Tisj 


KX LSJ requires the subroutines: LNORPC, LGAMPC. 
Declare variables. 


DOUBLE PRECISION DSEED 
REAL RPlsj, SDIisj, E#sj, Tisj, WN, 2B 


Calculate the Lambda Sigma Jump error fluctuation 
CElsj>). The error fluctuation is normally distributed 
with mean zero and standard deviation = SDlsj. 


CALL LNORPRC<(DSEED. Ne 
Elsj=SD1isjxN 


Calculate the time to the fluctuation (Tlsj). The time 
between flucuations is exponentially distributed with 
mean = 1/RPlisj. 


CALE EGAMNPCG (DSEBU ey tele 


Tlsj=E/RPlsj 
END 
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Wem emer emo eNO NCRGIe CC CCC CC CC CCC CC 
SUBROUTINE FALTR 

We@mtmemi ewe SC Cate CCC GC.C C.c. CC CCC CC 

SUBROUMGTNE  PALTRCDSEED sRPitr, LyeXf, Yf, TIf> 


Variables passed to FALTR are: DSEED, RPftr, L 
Variables returned by FALTR are: Xf, Yf, TTf 


* FALTR requires the subroutines: LRNDPC, LGAMPC 
Declare variables. Dimension arrays. 


DEVUBER ERESI SION DSEED 
Renew iitto Leenks, YF, Lif, U<2>, # 


Compute the x and y positions of the false transient 
Weis). Thee ecsition is uniformly distributed over 
the search area. 


CALL LRNDPCCDSEED, U, 2) 
XL=LAU C1) 
Y£f=LAU C2) 


Compute the time to the next false transient (TIf). 
Pie iie between ia)ce =r ranstents 1c ,exponentially 
fewer oUcetewieh mean = "l7KPit, in hours. 


CALE LOAM CDSEED, EE, 2. 1.0 


Lit -EARPL EE 
END 


oo 


C 
C 
C 
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CcCcccccchOCCCCCC CC CC COC Cee OC eee eee 
SUBROUTINE POSIT 
CccCcCccccccmcmcCCCCcUCCCCOC CCC CC CCC Cee 
SUBROUTINE POSIT(T, At, Yt, Ct, St, Xs, Ys, Csi 
Variables passed to POSIT are: 7T, At, Yt, Ct, Sli 
Ys, Cs, Ss 
Variables returned by POSIT are: At, Yt, Ks, Ys 
xX POSIT updates the variables Xt, Yt, Xs, Ys 
Declare variables. 
REAL T, AU, Ye, Ct, 2 St, ts 96s, CS eos moe 
Assign the value for Pl. 
PI=3. 1415926 


Compute target's x position (Xt) and y position jaa 


Xt=Xt+ (TKSt*SIN (CtxP1I/180)) 
Yt=Yt+ (TAStACOS (Ct *P1/180) ) 


Compute searcher’s x position (Xs) and y position 
CYS )% 


XS=Xst (TKSS*KSIN (CS¥PI/180)) 


Ys=Yst+ (T*ASS*COS (CS*P1/180)) 
END 
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Sg) ey Oe) 


OO OOo 


OOo O2es@ Oe aeG) 


ce < 


Ce 


Tememestee set eweweweee CCC CCC CUCUeCecCccCcCce 

SuUpROUtINE REPLCT | : 
Te meememeneme © eC CC CCG CC CC CC Cec CC ; 
SUE nROUt Neer EePECT CL, Ct, Styext, Yt, Cr, Trd 


Neampiables passed to REFLCT are: L, Ct, St, Xt, Yt 
Variables returned by REFLCT are: Cr, Tr 


Declare variables. 

Pam etise ob, ooe, YC, Gry, Ir, Piltemsts, Sty, Ty, Tx 
Assign the value of Pl. 

PI=3. 1415926 


Compute the x component of target speed (Stx). Compute 
the yeeomponent of target speed (Sty). 


Stx=StxSIN (Ct*xP1/180) 
Sty=StxcOs (Ct *P1/180) 


. Compute the time to reflection (Tr) and the reflected 


course (Cr) for the five cases: 


ie spectal cases = €t=0, 90, 180, 270, 300 
2. 0<CEC90 

S- GVOCEECLEO 

At eOCeE<27 0 

Sees Ose E360 


Special cases. 


PmeecGte EO, O) GOR. <Ct BQ: 360)) THEN 
ie Coster) 7S y. 
Cr=160 
Golo 10 

END IF 

Pe Cet. EO. 90) THEN 
tie U-AG 7 OoCx 
Cr=270 
20 1NG ys s08, 

END IF 

irecCr 2 EO. 180) THEN 
ie = 1 t/OLy 
Cr=0 
GOTEs10 

PND SE 

PreGame, 270) THEN 
te=—-AL/Stx 
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OE @O) 


COG) 


Cr=90 
GOTO; 10 
END” EF 


Case 1. 


IF GGCt. «GT a0) WAND @eeCCtE | LT. SOD wIE 
Tly= (Ean) Zocy 
Tx=(L- x7 St = 
Tr=MIN(CTy, Tx) 

LE lr EOS tye ney 
Cr=180-Ct 
BSE 
Cr=360-Ge 
END IF 
GOTO 10 
END IF 


Case cd. 


IFCCCt .GT. 90) VTAND. <Ct JLT. 260)) ire 
Ty=-Yt/Sty 
Tx= (L-Xto/7Stx 
Tr=MIN(CTy, Tx) 
If Clip weeo 2 iy) THEN 
Cra CO-Ct 
ELSE 
Cr=360-Ct 
END IF 
GOTO 10 
END IF 


Case 3. 


IFC CCt .GT. 180) 2AND. “CU 21 27ers 
Ty=-Ye7ety 
Tx=-At/Stx 
Tr=MINCTy, Tx) 

PE Cle ee EO ty RE 
Cr=540-Ct 
EBESE 
Cr=360-Ct 
ENDFIF 
GOTO 10 
END IF 


Case 4. 


IFC(Ct .GT. 270) .AND. “Cte bl. S600 ae 
y= Cla 6) 7 cay: 
Tx=-Xt7sStx 
Tr=MIN(CTy, Tx) 
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IF(Tr .EQ. Ty) THEN 


Cr—540-CE 
PESE 
€r=360-Ct 
BND IF 
GOTO 10 
END IF 


10 END 
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CCC C CCC CC CCC C CG C Cher ere eee 
SUBROUTINE TRACK : 

cCcCcCcccccmcmCOCCOCCCCCOCUCCCCCOCC CC € C CC Cee € Clee : 

SUBROUTINE TRACK(Xt, Yt, XKs, Ys, Ss, L, €s,maeee 


Variables passed to TRACK are: Xt, Yt, Xs, Ys, Ss, L 
Variables returned by TRACK are: Cs, TCs 


xX TRACK requires the subroutine REFLCT. 

Declare variables. 

REAL Xt, Yt, Xs, Ys, Ss, L, Cs, TCs, PlI, HOLD, Grae 
Assign the value for Pl. 


PI=3.1415926 


Gee Olé Compute the approach track course (Cs) -fomerae 
following cases: 
1. Special case where the y positions are eqiuem 
2 CCl CE wo mands. ( rt. . Gl) sys) 
3. (KES LEAKS) and eGyt 51. Ys) 
Aw CR bE NS) vance ey Gl ae 
De eure CE AS) (and “Cree -Li mes) 
-- Case 1. 


LECYSs EO. to wie 
IFCXt .GT. Xs) THEN 
Cs=90 
ELSE 
Cs=270 
END IF 
GOTO 10 
END IF 


Compute the common value for the angle of the track 
between the searcher and target (HOLD) to be useciie 
cases 2-0. 


HOLD= CATAN ( (Xt-Xs)/ CYt-Ys)))*180/7PI1 
Case e2. 


IF @QRt 2GEe XAs> - AND. (Yt AG Voor e esl 
Cs=HOLD 
GOTO 10 

END 
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Case 3. 

eecCxyt LE. Xs) «AND. CYT 
Cs=HOLD+180 
GOTe. 10 

END IF 

Case 4. 

PiGGCU wor. As) esAND. <Yt 
Cs=HOLD+360 
GOTO 10 

END IF 

Case 5. 

ieecKt «GE. Xs) . AND. CYt 
Cs=HOLD+180 

END IF 


wey 


Gl 


el. 


Call the REFLCT subroutine to 


course change (TCs). 


1s)) THEN 


Ys) ) THEN 


Ys) ) THEN 


Pomevne une co Next 


CAMP REPE@IGbh Ce,Se, Xs, Ys, Cr, Tr) 


TCs=Tr 
END 
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CCCCCCCCCC CC CC CC C CC Cee eee 
SUBROUTINE OUT ; 

CCCCCCCCCcCecCecCc CC CC CC Cee eee : 

SUBROUTINE OUIT<Na, Nz, Nm, Ftw, IW.) 1 max 


Variables passed to OUT are: Nd, Nz, Nm, Ftw, TW, 
Tmax, Td 
Variables returned by OUT are: none 


X OUT requires the subroutine HSCRN 
Declare variables. Dimension arrays. 


CHARACTERA12 FILE1 

INTEGER Nd, Nz, Nm, Ftw, ICOUNT 

REAL IW, Tmax, Td(€5000), Tdor(5000), SUNG suing 
SUMid, SUMim, SUM2, SUM3, SUM4, SUMS, SUM6, SUN’, 
XBAR, SD, SDx, XBARd, SBAR, RPe, RS, SUR(Ca00Ge 
PD (1 O0FFE Dt <1») 


Designate the output file. FILE1 —- holds the times to 
detection. 


FILE1=’ TIME. OUT’ 


Open FILE] for writing to. The other output 2 ame 
(SUBTRAN. DAT) is already open from the subroutine IN 
and is designated as device 20. 


OPEN (30, FILE=FILE]) 


Write to the user's screen a message indicating 
computing output statistics. 


WRITEC*, 14) 


Order the times to detection from smallest to largest 
(Terd). 


Tdor (1)=Td(1) 
DO 10 I=1,Nd 
DeeZ0e)=i-iatad 
IFCTaC(l)> SLI. Tder<)) THEN 
1der (i+ 1) =taorn)) 
1aor (3)—tarc) 
ELSE 
der (+1) =1a > 
GOTO 10 
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END IF 
CONTINUE 
CONTINUE 


Compute the mean (XBAR) and standard deviation (SD) of 
the times to detection. For the cases where detection 
did not occur before Tmax, assign a value of Tmax for 
the time to detection. Cases where the time to 
detection is zero are not included. Also include the 
905% confidence bounds on the estimate of the mean 
GSDs) 


DO 30 [=1,Nm 
Tdor (Nd+I)=Tmax 
CONTINUE 
SUMd=0 
SUMm=0 
SUM1d=0 
SUM1m=0 
DO 40 [=1,Nd 
SUMd=SUMd+Tdor (I) 
SUM1d=SUMid+ (Tdor (1) *xe) 
CONTINUE 
DO 50 [=1, Nm 
SUMm=SUMm+Tdor (Nd+I1) 
SUM1im=SUMim+ (Tdor (Ndt+I)xkx2) 
CONTINUE 
XBAR= (SUMd+SUMm / (Nd+Nm) 
SD=SQRT (( (SUM1d+SUM1m — ( (Nd+Nm) * (XBARK*2)))/(Nd+Nm-1)) 
SDx=1.645xXSD/ ( (Nd+Nm) kx. 5) 


Construct a empirical survival function (SUR) from the 
ordered time to detection data. Under the assumption 
that the data is from an exponential distribution, 
perform a linear regression on the natural log of SUR 
with respect to the ordered times Tdor. The slope of 
the resultant line is an estimate for the rate 
parameter of the exponential distribution (RPe). 1/kRPe 
is another estimate for the mean time to detection. 
This estimate is independent of the detections that 
occurred at time zero thus those values are not 
included. Also computed is the R-squared value of the 
linear regression. This value is an indication of how 
good the exponential assumption is. 


SUM2=0 

SUM3=0 

SUM4=0 

ICOUNT=0 

DO 60 I=1,Nd 
SUR (I)=1- (REAL (1)/7REAL (Ndt+Nm ) 
Peecoun Clo, LE. 0) THEN 
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60 


80 
70 


1 


90 


ICOUNT=ICOUNT+1 
SUMd=SUMd-Tdor (1) 
SUM1Id=SUMId- (Tdor (1) **2) 
GOT@sso 
END Re 
SUM2=SUM2+ALO0G (SUR CI) ) 
SUM3=SUM3+ (Tdor (1) XxALOG (SUR C1))) 
SUM4=SUM4+ CALOG (SUR C1) ) **2) 
CONTINUE 
DONO t= 1 ae 
DO 80 J=1,Nd 
TF<C1-SURC)> (GE. (REAL CID 710 ae 
PDG1)=1=sURGI> 
FOE Ciw—Tdor cs 
CO1@s 70 
END IF 
CONTINUE 
CONTINUE 
XBARG=SUMd/ (Nd- I COUNT) 
SBAR=SUM2/ (Nd-ICOUNT) 
RPe= (SUM3— (XBARGAXSUM2))7 
(SUM1id- ((Nd- ICOUNT) * (CXBARG*k2))) 
SUMS5S=0 
DO 90 I=1,Nd 
IF(SUR(I)D . LES O) GOlGs se 
SUMS=SUMD+ ( (Tdor t)) -XABARd) * (ALOG (SUR (C1) =SBAR ye 
CONTINUE 
SUMO=SUM1d-— ((Nd- ICOUNT) * (XBARd* Xe) ) 
SUM7=SUM4—- ( (Nd- ICOUNT) * (SBARKX2) ) 
RS= (SUM54*2) 7 (SUM6XSUM7) 


Write the output to the screen and the two output data 
files. 


WRITE(*,114) 
WRITE (*, 24) 
WRITE(20, 24) 
WRITE(*, 34) Nz 
WRITE(20, 34) Nz 
WRITE(*, 44) Nm 
WRITE(20, 44) Nm 
IF(Ftw .EQ. 0)THEN 
WRITE(*,54)TW 
WRITE(20,54)TW 
ELSE 
WRITE (*, *) 
END IF 
WRITE (*, 64) XBAR, SDx, SD, (-1/RPe), RS 
WRITE(20, 64) XBAR, SDx, SD, (-1/RPe),RS 
WRITE (*,74) 
WRITE (20,74) 
DO 100 1=1,5 


108 


Vii ea Clr bl <)>) POC lL+5), PDI (1+5) 
Wheelie (20.34) Pore werd Gi FDC]+5), PDTCI+5) 
100 CONTINUE 
Wiel DE CX. x*) 
CALL HSCRN 
WRITEC(X, 94) 
DO 110 [=1, (Ndt+tNm 
WRITE CS0, *>) Tdor <i) 
LtO CONTINUE 
WRIdtE Cx, 104) 


CALL HSCRN 

C 

Cae Write format statements. 

C 

14 ihe czs (7 71k, COMpuLing Output statistics... .'/) 

24 FORMA Czex, 26 t ')726X,°* SUBTRAN OUTPUT DATA *’/ 
i Ee Gar | gag al, GE 

34 PonnAI Cex NO. OF REPS. WITH@ DETECTION TIME OF ', 
i TNO ee a ey LAD 

44 PORMA) (Oye eeeer REPS. WETH DETECTION TIME’, 
1 pCOkEiwger THAN Tmax ..8.0%. ry aD 

54 eke lA, SALECULATED TRACK WIDTH VALUE ', 
it Os RE OI ee (5 18 ae, 

64 MORMAT<7/2ze., TIME TO DETECTION STATISTICS’ /26X, 
1 AO = wee rcOn, Mean Se eee YS, 
1 poe So feed 720%, ota. Dev, = ',F%.2/ 
1 Zee eeabambdawe—  §,F/’.2/20%,  R-Squared = ', 
1 eae) 

74 Pee CeO emieouweat Det.) le/7,° By time’, T47, 
i SrTihemer moet. 104, By time /110,13¢'-"), 
1 Wey een 1a LOO =), 164, 7 C=") 

84 POW ne Mi acme 2c, Oo. 1, 152, F3.1, 165,F6.1) 

94 Poke ie, Waiting CO Output file....° 7?) 


104 FORMAT (4X, ’Output data is located in the following ’, 
Peis a vee oD TRAN DAL = input *, 
"parameters, time to detection '’, 
Petal Sule owme, (sy, 7.0 LIME GUT aelLLmeSs eto. a, 
meevecLions..477) 

114 FORMAT (23 (7) ) 


END 


a 
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CCCCCCCC C&C CC CC C © Cee ene eee eee ee 
SUBROUTINE HSCRN 

CCCCCCcCcCe CCC CCC CC C Chee ee ee eee 

SUBROUTINE HSCRN 


Write to the user's terminal a "how to continue” 
message. 


WRITE CX, 14) 


Perform a generic READ statement to hold the screen, 


READ Cx, X) 
Write format statement. 


FORMAT (26X,'- Press <ENTER> to continue —-') 
END 
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_Zz=7—OMeeeee oes eG SC CC CCC Ce CCC CCCCCCG 
C 
C SUBROUTINE LRNEPG 
C 
femme Gc © GCeC CC ee CC CCcCcCCOCcCcmmecmmcCrmlCcCC CCC Ce 
C 
SUBROUTINE LRNDPCCDSEED, U, WN) 
C 
INTEGER*4 N, I 
REAL UCN) 
DOUBLE PRECISION D3iM1, DSEED, D3i 
C 
C a3imil = 2xx31-1 
C a31 = 24%*31 
C 
DATA D31M1/72147483647. DO/Z 
DATA D31/2147483648. DO/ 
C 
DO 5 I=1,N 
DSEED = DMOD(16807.DOXDSEED, D31M1) 
DiGi ee USEED,D31 
5 CONTINUE 


END 
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ccccccccccCceCce Cc CC C € € CC Cee ere 
C C 
C SUBROUTINE LNORPC C 
C | C 
Cccccccceccececeecce Cc ce € C € Cle Clee. 
C 

SUBROUTINE LNORPCCISEED, A, WN) 


This program will generate a vector of normal random 
variables using according to the sine-cosine method. 


Oa 


DPNEEGER. AN, PN 

REAL ACN) ,.U<(2)7'S5, Ww; UL, 2.SesR 
DOUBLE PRECISION “[SEED, Et 
DATA PI/3.14159265358979D0/ 


C DATA IND/1/ 
IND = 1 
DO 100 I=1,N 
IND = -IND 
IF CIND. GEO oe Geter 
10 ‘GALL LRENDPC GISEED MUR 27 
© =SQRT (-2xALOG CU C1))) 
W = 2*¥PTAU C2) 
ASTAR = S*COS (CW) 
ACI) = SKkSINCW) 
GOT@ 109 
ZO ACI) = XSTAR 
100 CONTINUE 
RETURN 
END 


BY 


C 
C 
@ 
C 
C 
C 


cS ee omen ome 
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C 
o 

SUBROUTINE LGAMPC E 
| e 

C 


Cemememememerermere CeCe ©€ CCC CeCe CCcCccccce 


SUBROUTINE LGAMPCC(DSEED, G, N, A) 
DEUBLCE PRECISION DSEED 

DIMENSION G(N), UN(C2) 

DAIWA ASAVE 7-1.7,XLL/-T1.7 

Peewee tee? CO TO 1500 

GAMMA GENERATOR A>1. SCHMEIZER’S ALGORITHM 
IF (A .EQ. ASAVE) GO TO 50 

ASAVE = A 

X1=0 

X2=0 

F1=0 

F2=0 

Se Aas «1. 

RSX=X3-X3XALOG (X3) 

D = SQRT (X3) 

Pee Goren) Ge TO 10 


M2 ko, = D 

nae — ones Cl. — 1. 7D) 

Rieke => 1. = AS/XKI 

Fil = EXPCXSKXALOGC(X1) —-X1+R3X) 
Fo = EXPCXS*ALOG(X2) -— Xe + RSX) 
X4 = X3 + D 

X5 = X4*C1. + 1./7D) 

Moke = le “= ekS7XSD 

F4 = EXP CX3*ALOG(X4) - K4 + RSX) 
Row EXE CXS tALOG CXS) =. AD + RSX) 
Pie = F2kOKS -hoD 

Pe = F4X* (X4-X3) + Pil 

ie = Pit (X2-K1) > 4 Pe 

P4 = FOS* (X5-X4) + PS 

ie, = AS=AZ —-P1l +). P4 

RO = A4-AS-P2tP1 + PS 

Momo ec PloA CXS KIDD + PS 

Ro = (F4—-Fo ACAD k4)0%.5 + P7 

Ro = =FPIVxXLL + Fé 


PaOe-shS7 7k + FS 

ireneecOr Ar molneee., 14, .otilebe F4 b>, PP], Pe, P3,P4,P5,P6, 
East or LO 

ROKMAT (SF20.6) 

DO 1400 [{!f=1,N 

CALL URNDPC(DSEED, UZ 1) 

U = UXP10 

heals GisP4) GO TO 500 

Peeve Gi.arid GO FO 200 

X = Xe + U/Fe 

GO TO 1400 
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200 


300 


400 


500 


600 


700 


800 


900 


1000 


1300 
1400 


eno 
LSOG 


BUSS, 


IF (U.GT.P2) GO TO 300 

X = X38 + (U-P1)/F4 

GO 1@) 2400 

IF (U.GT.P3) GO TO 400 

X = X1 + (U-P2)/F1 

GO TO 1400 

X = X4 + (U-P3)/F5 

GO TO 1400 

CALL LRNDPC(DSEED, W, 1) 

IF (U.GT.P5) GO TO 600 

X S02 ites arena 

IF ((U-P4)/(P5-P4) .LE. W) GO TO 1400 
V=F2 + (U-P4) / (X3-X2) 

GO TO 1300 

IF Ui, GivPe, (co slemaae 

X=X3+ (X4-X3) *W 

IF ((P6-U)/(P6-P5) .GE. W) GO TO 1400 
V = F4 + (U-P5)/(X4-X3) 

GO TO 1300 

IF (U .GT. P8&) GO TO 900 

CALL LRNDPC(DSEED, W2, 1) 

IF (W2 .GT. W) W=W2 

IF (U.GT.P7) GO TO 800 

X=X1+ (X2-X1)*W 
V=F1+2. kWk (U-P6) / (X2-X1) 

IF (V.LE.F2*W) GO TO 1400 

GO TO 1300 

X=X5 - Wk(X5-X4) 

V = F5+2. kWk (U-P7) / (X5-X4) 

IF (V.LE.F4*W) GO TO 1400 

GO TO 1300 

IF (U.GT.P9) GO TO 1000 

U= (PO-U) / (P9-P8) 

X=X1-ALOG (U) /XLL 

IF (%.LE.0) GO TO 260 

IF (W.LT. (XLL*(X1-X)+1.9/U) GO TO 1400 
V = WKF1XU 

GO To) 1300) 

U = (P10-U)/(P10-P9) 
X=X5-ALOG (U) /XLR 

IF (W.LT. (XLRK(X5-X) + 1.)97U) GO TO 1400 
V = WKF5xU 

IF (ALOG(V) .GT. X3*ALOG(X) + R3X - X) GO TO 100 
EO Ye 

RETURN 

GAMMA GENERATOR A<=1. ALGORITHM GS. 
CONTINUE 

IF (A .EQ. ASAVE) GO TO 1550 

ASAVE = A 

AINV=1/A 

B=A*.36787944 + 1. 

CONTINUE 
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1560 


1600 


e700 


1800 


1810 
1900 


Poon to Ge. 10) LeO0 
DGC OO l=1,N 

GOnTInNuUE 

GAPE SER NDEC CDSEED, UN, 2? 
EXPL=—-ALOG (UN (1) ) 

P=BXUN (2) 

eC Otel oe GO, TO 1600 
X=PAXKAINV 

Die. Glo eexPi) GO TO S60 
G Gl) -—x 

Cemt0 1700 

CONTINUE 

X=-ALOG (CB-P) XAINV) 

Wee CIA FALGG (AD .GT. EXPL) «GO TO 1560 
GCID=X 

CONTINUE 

GOTO 1900 

CONTINUE 

CAEL SMERNOreCDSeED, G, WN) 
DO 1810 I!I=1,N 

Clo eae Rh LOG CG Cl) 
CONTINUE 

CONTINUE 

RE PUR 

END 
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